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Abstract 

The U(l) extended supersymmetric standard model (UMSSM) can accommo¬ 
date a Higgs boson at 125 GeV without relying on large corrections from the 
top/stop sector. After imposing LHC results on the Higgs sector, on H-physics 
and on new particle searches as well as dark matter constraints, we show that this 
model offers two viable dark matter candidates, the right-handed (RH) sneutrino 
or the neutralino. Limits on supersymmetric partners from LHC simplified model 
searches are imposed using SModelS and allow for light squarks and gluinos. More¬ 
over the upper limit on the relic abundance often favours scenarios with long-lived 
particles. Searches for a Z' at the LHC remain the most unambiguous probes of 
this model. Interestingly, the D-teim contributions to the sfermion masses allow 
to explain the anomalous magnetic moment of the muon in specific corners of the 
parameter space with light smuons or left-handed (LH) sneutrinos. We finally em¬ 
phasize the interplay between direct searches for dark matter and LHC simplified 
model searches. 
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1 Introduction 


The discovery by the ATLAS and CMS collaborations [1-3] of a 125 GeV Higgs boson 
whose properties are compatible with the standard model (SM) predictions conpled with 
the fruitless searches for new particles at Run I of the LHC [4-7] has left the community 
with little guidance for which direction to search for new physics at the TeV scale. The 
dark matter (DM) problem remains a strong motivation for considering extensions of the 
SM, in particular supersymmetry. 

In the minimal supersymmetric standard model (MSSM), the Higgs couplings are to a 
large extent SM-like, especially when the mass scales of the second Higgs doublet and/or 
of other new particles that enter the loop-induced Higgs couplings are well above the 
electroweak scale. This is to be expected in any model where the Higgs is responsible for 
electroweak symmetry breaking. The main challenge for the MSSM is however to explain a 
Higgs mass of 125 GeV. To achieve such a high mass requires large contributions from one- 
loop diagrams involving top squarks — in fact the loop contribution has to be of the same 
order as the tree-level contribution — thus introducing a large amount of fine-tuning [8,9]. 
The fine-tuning is reduced in extensions of the minimal model containing an additional 
singlet scalar field [10-13]. For example in the next-to-minimal supersymmetric standard 
model (NMSSM), terms in the superpotential give an extra tree-level contribution to 
the light Higgs mass, thus reducing the amount of fine-tuning required to reach ruh = 
125 GeV. A doublet-singlet mixing can also modify significantly the tree-level couplings of 
the light Higgs. In the UMSSM, where the gauge group contains an extra U(l) symmetry, 
contributions from U(l) D-terms in addition to those from the superpotential present in 
the NMSSM, can further increase the light Higgs mass [14,15] reaching easily 125 GeV 
without a very large contribution from the stop sector. Furthermore, because the singlet 
mass is driven by the mass of the new gauge boson which is strongly constrained by 
LHG searches to be above the TeV scale [16,17]\ the tree-level couplings of the light 
Higgs are expected to be SM-like, in agreement with the latest results of ATLAS and 
GMS [19,20]. This heavy Z' was also found to increase the fine-tuning of supersymmetric 
models with U(l) extended gauge symmetry [ 21 ]. Another nice feature of the UMSSM 
(as the NMSSM) is that the p parameter, generated from the vacuum expectation value 
(vev) of the singlet field responsible for the breaking of the U(l) symmetry, is naturally 
at the weak scale. Finally, this model is well motivated within the context of superstring 
models [22-26] and grand unified theories [27,28] . 

The range of masses for the Higgs scalars and pseudoscalars were examined in a variety 
of singlet extension of the MSSM [15]. The parameter space of a similar model with a 
new U(l) symmetry, the constrained EqSSM, compatible with the Higgs at 125 GeV as 
well as limits on the Higgs sector and providing a dark matter candidate was examined 
in [29]. In this model the RH sneutrino does not have a U(l) charge and is expected to 
be very massive. The h ^ 77 signal in a U(l) extended MSSM model was discussed 
in [30,31] with emphasis on the region that leads to an increase in the two-photon signal. 
Additional non-standard decays of Higgs particles were found in [32,33]. However, since 
the mass of the additional singlet Higgs is expected to be very large due to strong limits 
on the Z' boson mass, it does not affect the property of the lightest Higgs which is hence 
expected to be SM-like. In the UMSSM model considered here RH sneutrinos can be 
charged under the additional U(l) symmetry, hence this model gives a new viable dark 

^In this paper we concentrate on a Z' above the electroweak scale, for scenarios with light Z' see [18]. 
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matter candidate in addition to the lightest neutralino as observed in [34] . The properties 
of a RH sneutrino DM were also examined in the U(1)b-l [35] and U{1)b-l x U{1)r 
extensions of the MSSM [36] . Note that in such models the sneutrino vev’s were found to 
play an important role in the vacuum stability [37]. Furthermore the Z' can contribute 
to the stabilization of the Higgs potential [38] . 

In this paper we explore the parameter space of the UMSSM (derived from E^) that 
is compatible with both collider and dark matter observables. We include in particular 
the Higgs mass and signal strengths in all channels, LHC constraints on Z' and on super- 
symmetric particles, new results from H-physics, as well as the relic density and direct 
detection of dark matter. Specifically we take into account the most recent LHC results 
for supersymmetric particle searches based on simplified models using SModelS [39,40]. 
This allows us to also highlight the signatures not well constrained by current searches 
despite a spectrum well below the TeV scale. One salient feature of the model is that 
large D-term contributions can significantly reduce the mass of RH squarks thus splitting 
the u-type and d-type squarks and weakening the constraints on first generation squarks. 
Another feature, which is also found in the MSSM, is that the relic density upper limit 
favors a neutralino with a large higgsino or wino component as the lightest supersymmet¬ 
ric particle (LSP). Scenarios with a higgsino LSP can easily escape current search limits. 
For example simplified model limits from top squark searches rely on the assumption that 
one decay channel is dominant, while for higgsino LSP branching ratios into and hxl 
can both be large, thus the mixed channels where each stop decay into a different final 
state are important. Since a higgsino or wino LSP may be associated with a chargino 
which is stable at the collider scale, we also impose the DO and ATLAS limits originating 
from searches for long-lived particles. On the remaining parameter space, we then discuss 
the expected spectra of SUSY particles, the expectations for the signal strengths for the 
Higgses as well as dark matter observables in direct and indirect detection. 

In general we do not attempt to explain the observed discrepancy with the standard 
model expectations in the muon anomalous magnetic moment. However, we highlight 
the region where the model can explain this discrepancy and investigate how it may 
escape simplified model limits from the LHC. The interplay between the muon anomalous 
magnetic moment constraint, LHC and DM limits was recently studied in the MSSM [41]. 

In contrast to previous studies [34,42] we explore the impact of LHCSTeV results on 
Higgs and new particle searches from the 8 TeV run on scenarios with arbitrary U(l) 
originating from Eq. Moreover we consider both the cases of a neutralino and a RH 
sneutrino dark matter. We further examine the implications of dark matter searches in 
these scenarios. An attractive feature of the model is the possibility to obtain ruh = 
125 GeV despite small values of tan/5. The phenomenology of Higgs and SUSY searches 
could thus differ from that of the much-studied MSSM. 

This paper is organised as follows. In section 2 we describe the model. Section 3 is 
devoted to a detailed view of the Higgs sector. Section 4 presents the different constraints 
used in our study. Section 5 contains the results for several sectors of the model after 
applying a basic set of constraints mostly related to Higgs and R-physics observables and 
after applying the DM relic abundance limits. Section 6 is dedicated to the application 
of the LHC simplified models searches on the remaining allowed parameter space of the 
UMSSM. Section 7 contains a summary of LHC constraints after Run I and suggestions 
on how to extend simplified models searches to further probe the model. Section 8 shows 
prospects for probing the Higgs sector and section 9 prospects from astroparticle searches. 
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Our conclusions are presented in section 10. 


2 The model 

The symmetry group of the model is SU{3)c ® SU{2)l ® U{1)y ® U{iy and we assume 
that this model is derived from an underlying Eq model. In this case the ?7(1)' charges of 
each field F of the model are parameterized by an angle 6eq as 

Q'p = cos Ose Q!^ + sin Q'^, (2.1) 

where 6eq G [—7r/2,7r/2] and the charges and Q'^ are given in Table 1 for all fermionic 
fields that we will consider [43,44]. The dependence on 6eq of the 17(1)' charge of some 
matter fields is shown in figure 1. 

The matter sector of the Eq model contains, in addition to the chiral supermultiplets of 
the SM fermions, three families of new particles, each family containing : a RH neutrino, 
two Higgs doublets (i7„,i7rf), a singlet, and a colour SU{3)c (anti)triplet. While the 
complete matter sector is needed for anomaly cancellations, for simplicity we will assume 
that all exotic fields, with the exception of three RH neutrinos, two Higgs doublets and 
one singlet, are above a few TeV’s and can be neglected. Similarly in addition to the 
MSSM chiral multiplets we will only consider the chiral multiplets corresponding to these 
fields, that is the multiplet with a singlet S and the singlino S and another multiplet with 
RH neutrinos UiR {i G {e,/r, r}) and their supersymmetric partners, the sneutrinos, DiR. 


Q'q Q; Q:, Q'e Q' Q'e Qk Qk gg 

-1-13 3 -5 -12 -2 0 

111111-2-24 

Table 1: ?7(1)' charges of all matter fields considered. 



Figure 1: ?7(1)' charges of some matter fields in the UMSSM as a function of 
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Finally the UMSSM model contains a new vector multiplet, with a new boson B' and 
the corresponding gaugino 5'. The superpotential is the same as in the MSSM with /i = 0 
but has additional terms involving the singlet, 

WuMSSM = VVMSSM|/i=o + ^SHuHd + + 0(TeVs) (2.2) 

where is the neutrino Yukawa matrix. The vev of S, (S) = ^ breaks the U{iy 
symmetry and induces a jjl term 

/r = A^. (2.3) 

Note that for Oeq = 0 the U{iy symmetry cannot be broken by the singlet field since 
Q'g = 0. Note also that the invariance of the superpotential under 17(1)' imposes a 
condition on the Higgs sector, namely + Q's — 0- The soft-breaking Lagrangian 

of the UMSSM is 

“^UMSSM ='=^MsWlfe=0 ~ T + h.C.^ — 

- ml\sy - (XAxSHuHd + h.c.) + (!)(TeVs), 

with the trilinear coupling A\, the B' mass term M[, the singlet mass term ms- The 
soft sneutrino mass term matrices and are taken to be diagonal in the family 
space. Note that our study is based on the UMSSM model with parameters defined 
at the electroweak scale, we make no attempt to check the validity of the model at a 
high scale. We now describe briefly the sectors of the model that will play a role in the 
considered observables. 

2.1 Gauge bosons 

The two neutral massive gauge bosons, and Z' = B' can mix both through mass 
and kinetic mixing [42,45]. In the following we will neglect the kinetic mixing^. The 
electroweak and U(l)' symmetries are broken respectively by the vev’s of the doublets, 
u„/-\/2 = {Hu), = {Hd} and singlet, Vsl\f2 = {S). The mass matrix reads 


f M|o A| \ 

V A| M|, J ’ 


(2.5) 


where 


= J^(Vn + ^d) 

Ml = + + 


CA = - q'„A) 


2o 


W 


( 2 . 6 ) 

(2.7) 


^The impact of the kinetic mixing on the Higgs boson mass and on the Z' and DM phenomenology 
was examined in the U{1)b-l extension of the MSSM in [35,46,47]. 
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where g 2 = ejsw, Qi = \/5/3pi, gi = e/cw and cw (syi/) is the cosinus (sinus) of the 
Weinberg angle. Diagonalisation of the mass matrix leads to two eigenstates 

Zi = cos azZ^ + smazZ' 

Z 2 = — sm.azZ^ -\-cosazZ' ( 2 . 8 ) 

where the mixing angle is dehned as 


sin 2az 


2A| 
U\ - 


(2.9) 


and the masses of the physical helds are 

+Mf, T ■ (2.10) 

Precision measurements at the Z°-pole and from low energy neutral currents provide 
stringent constraints on the Z^ — Z' mixing angle. Depending on the model parameters 
the constraints are below a few 10“^ [48,49]. The new gauge boson Z 2 will therefore have 
approximately the same properties as the Z'. As input parameters we choose the physical 
masses, Mz^ = 91.187 GeV, Mz 2 and the mixing angle, az- From these together with the 
coupling constants, we extract both the value of tan/3 = v^/vd and the value of Vg. Note 
that as in [50] we adopt the convention where both A and tan (5 are positive while /r (and 
then Vs) and A\ can have both signs. From eqs. (2.7) and (2.9), 


cos^ /3 


Q\ 




Q'hs 


( sin2Q;z(M|^ - M^Jcw 

V v‘^9'i92 



( 2 . 11 ) 


where + nj. 

For each 17(1)' model the value of tan/3 can be strongly constrained as a consequence 
of the requirement 0 < cos^/3 < 1. For example for the U{1)^ case with sina^ > 0 and 
Mz 2 3> Mzi the value of tan/3 has to be below 1. The reason is that for this choice of 
6*£g we have 

^2 ^ 929[ 

C\Y \/24l 

For other choices of parameters the value of tan/3 can be very large, 0(100). Another 
interesting relation is found for the case of small mass mixing between and Z' namely 

^ In this limit /3 is determined from the ?7(1)' charges only, 

^2 


(tan^ (3 — l)'u^ < 0. 


( 2 . 12 ) 


cos^ /3 ~ 


Q) 


Hu 


Q'h. + Q) 


Hu 


(2.13) 


One might think that small values of tan /3 are problematic for the Higgs boson mass 
since the MSSM-type tree-level contribution becomes very small. However, as we will see 
below, additional terms to the light Higgs mass and especially their dependence on az 
can help raise its value to 125 GeV. 
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2.2 Sfermions 


The important new feature in the sfermion sector is that the U{iy symmetry induces new 
D-term contributions to the sfermion masses. These are added to the diagonal part of 
the usual MSSM sfermion matrix, and read 


= yiQ!F (Q'i 




Q', 




+ 




(2.14) 


where F G {Q, u, d, L, e, f}. 

For large values of Vg the new T>-term contribution can completely dominate the 
sfermion mass. Moreover this term can induce negative mass corrections, even driving the 
charged sfermion to be the LSP. Thus the requirement that the LSP be neutral (either 
the lightest neutralino or RH sneutrino) constrains the values of Oeq (unless one allows 
large soft masses for the sfermions). For example, for — tan“^(3A/3/5) < Oeq < 0, the 
corrections to the d-squark and to LH slepton masses are negative, while for 0 < Oeq < 
tan“^(-y/3/5) the corrections to the u-squark and RH slepton masses are negative. The 
latter implies that the u-type squarks (and in particular the lightest top squark) and the 
RH sleptons can be the Next-to-LSP (NLSP). Interestingly for Oeq = — tan“^(3y/3/5) ~ 
— 1.16 the LH smuon/sneutrino can be sufficiently light to contribute signfficantly to 
the the anomalous magnetic moment of the muon and bring it in agreement with the 
data [51,52]. 


2.3 Neutralinos 

In the UMSSM the neutralino mass matrix in the basis (R, IT^, Ha, Hu, S, B') reads {c^ = 




= sin P) 

Ml 

0 

— MzoCpSw 

MzoSjjSw 

0 

0 

0 

M 2 

MzoCjsCw 

— MzoSpCw 

0 

0 

— MzoCjjSw 

MzoC/^Cw 

0 

-n 

\ Vu 

^V2 

Q'Hd9l^d 

MzoS^Sw 

— MzoSjsCw 


0 

^V2 

Q'huSi'^u 

0 

0 

_ 

\ Vd 

0 

Q's9[^s 

0 

0 


Q'hu9i'^'k- 

Q'sdl'^s 

M[ 


\ 


M^o = 


(2.15) 

Diagonalisation by a 6x6 unitary matrix leads to the neutralino mass eigenstates : 

e {1, 2,3,4, 5, 6}. (2.16) 

The chargino sector is identical to that of the MSSM. 

Several studies have analysed the properties of the neutralino sector in the UMSSM 
[53,54], in particular as concerns the neutralino LSP as a viable DM candidate [55,56]. In 
the weak scale model, the LSP can be any combination of bino/Higgsino/wino/singlino 
and bino’. However, as we will show, the LSP is never pure bino’, the pure bino and 
singlino tend to be overabundant while pure higgsino and wino lead to under abundance 
of DM. 
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3 The Higgs sector 

The Higgs sector of the UMSSM consists of three CP-even Higgs bosons hi,i e {1,2,3}, 
two charged Higgs bosons and one CP-odd Higgs boson A^. 

The Higgs potential is a sum of F-, D- and soft supersymmetry breaking-terms be¬ 
longing to the UMSSM Lagrangian : = Uf -h Vd -l- Uoft) where 


Uf = \XH^ • + |A5|2 {\Hd\^ + \H^\^) , 

Vd = 




9? 


{\Hd\^ - \H^\y + ^ {\Ha\^\Hu\^ - \Hu • Ha\^) 

+ Y + Q!H^Hu? + ^s\S?f , 

Usoft = m]j^\Hd\^ + + ml\S\^ + {A^XSH^ • + h.c.). 

At the minimum of the potential the neutral Higgs fields are expanded 


(3.1) 


as 


1 


^d = i^d + <Pd + Wd) , H^ = + (t)u + Wu) , 5' = ^ (^3 -F cr -F zO , (3.2) 


^ V^a -ra ^ ^ 

while the charged Higgs : 

Ha = — cos l3Gw- + sin I3H , = sin f3Gw+ + cos ^ 

with Gw the Goldstone boson. 

The minimization conditions of l{ygg are [15] 


(3.3) 


«•)' = - 


9i + 92 


0'2 fP 
^h,9i 


1 


9i + 92 


,/2 




[A^ -F QHaQs9i] + 


XAxVsVy 

V2vd 


KD = 2 


’’fi'l + 92 


-X^-Q'hAJi 


Vd-^ 


9i + 92 


+ QU? 




1 


1 


XAxVuVd 


— - 2 + Q^HiQ's9i] '^d~ 2 Q'HnQ's9i] 3 -^ 


(3.4) 


The tree-level mass-squared matrices for the CP-even (A4°) and CP-odd (A4°) Higgs 
bosons can be written in the basis {i?°, i?°, S} using the relations 


(A4°) 




d4'id<l)j 


0 




difidifj 


0 


(3.6) 


where {4>\,<p2A^) = (<^a.<(*«><t) and (<pi,<P 2 ,“Pa) = (ted.ttu,!;). For the neutral CP-even 






















Higgs bosons the relations are 


(•«"+)„ = 


9i + gl 


+ ^kg'i 




4 

Qi + ol 




XAxVsVy 

V^Vd 


- - Q'hAJi 


VuVd - 


XAxVs 

V2 


{^1)iS=[>^" + Qh,Qs 9?] VsVd 


(Ml),, = 


gl + gl 


+ Q'ka'i 


(■^ + )23 “ + ^Hu^Sg'k 

XAxVuVd 


XAxVu 

V2 

XAxVsVd 

\/2vu 

XAxVd 




For the CP-odd sector the mass matrix 


= 


XAx 

75 


.V2 


( VsVu 
Vd 


Vs Vu 


\ 


VsVd 


Vu Vd 


Vd 


VuVd 


leads to 


«A = 


XAxV2 f ■ 2 na 


The charged Higgs mass at tree-level reads 


-|- 


XAxV2 X^ 


sin 2/3 


(3.6) 


(3.7) 


(3.8) 


(3.9) 


The radiative corrections to the Higgs sector are given in appendix A. 

The lightest Higgs is usually SM like but can be heavier than in the MSSM. Indeed 
the tree-level lightest Higgs boson mass squared, which can be approximated by [57] 


^hi , tree 


~ M|o cos^ 2/3 -|- ^A^u^ sin^ 
AS2 / Aa sin^ 2/3 

“ V 7 


2/3 + pf [Q'h^ cos^ P + Q!h^ sin^ 

+ Q'hu /^) , 


(3.10) 


receives three types of additional contributions as compared to the MSSM. The first one 
proportional to A is also found in the NMSSM, the second one comes from the additional 
U(l) gauge coupling g[ and the last arises from a combination of pure UMSSM and 
NMSSM terms. The first term is not expected to play as important a role as in the 
NMSSM since A is small. This is because A is inversely proportional to the vev of the 
singlet Higgs field which is in turn related to the mass of new gauge boson, see eqs. (2.3) 
and (2.6). The strong dependence of the latter two terms on the U (1)' charges means that 
the size of the tree-level contribution to the Higgs mass will mostly depend on the value 
of Oeq ■ We illustrate this taking the limit of small mass mixing between the two Z bosons 
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as given in eq. (2.13). Figure 2 shows that in this limit does not exceed the MSSM 
upper bound and that its value depends strongly on 9eq- The maximum is reached for 
|6*Eg| = tan“^(-\/3/5) 0.66 which corresponds to cos^/5 or sin^/5 = 1. A smaller value 

for the maximum m^’’®® is found for lower values of Mz 2 since the last term in eq. (3.10) 
then gives a larger negative contribution to the tree-level mass. Furthermore, the tree- 
level mass tends to be suppressed for small fi values since these are linked to small values 
of A and thus to a small contribution from the NMSSM term. This behaviour shown in 
figure 2 is mostly observed for \fi\ < 2 TeV and |Aa| < 4 TeV. In the general case with 
non-zero — Z' mixing, a large tree-level contribution can be obtained for a wider range 
of parameters and a mass of 125 GeV for the lightest Higgs boson can easily be reached 
even with a small contribution from one-loop corrections that comes predominantly from 
the stop/top sector, as in the MSSM. 



2 

Figure 2: The tree-level light Higgs mass in the approximation olz <C as a function 

^2 

of 6eq for different values of |Aa|, Mz 2 and |/i|. 

Typically the Higgs spectrum will consist of a standard model like light Higgs, a 
heavy mostly doublet scalar which is almost degenerate with the pseudoscalar and the 
charged Higgs, and a predominantly singlet Higgs boson. The latter can be either h 2 or 
ha, depending on the values of the free parameters of the model, in particular Mz 2 and 
A\. The singlet Higgs is never hi because its mass depends on Vg which is large due to 
the lower bound on Mz 2 - 


4 Constraints on the model 

4.1 Higgs physics 

For the Higgs sector we require that the light^ Higgs mass lies in the range ruh^ = 

125.1 ± 3 GeV allowing for a theoretical uncertainty around 2 GeV. We impose con¬ 
straints on the Higgs sector keeping only points allowed by HiggsBounds-4.1.3 [58] and 
by HiggsSignals-1.2.0 [59] at 95% G.L (p-value above 0.05). We also use constraints 

^ Strictly speaking it is also possible that the Higgs at 125 GeV corresponds to h 2 , however we did 
not find such points in the scan. 
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contained in NMSSMTools [60], in particular the one on the heavy Higgs search in the 
decay mode that rules out some of the large tan (5 region. 

Note that the Yukawa couplings evaluated at the SUSY scale which enter the compu¬ 
tation of the Higgs boson masses must remain perturbative. We require that all Yukawa 
couplings stay below \/4^ at the SUSY scale. This condition will impose restrictions on 
both the small and the very large tan (5 values (recall that tan P is not a free parameter 
of the model). Yukawa couplings within the perturbative limit can nevertheless induce a 
very large width for some of the Higgs states, since we work in the context of elementary 
Higgs particles we impose the condition T{hi)/mh^ < 1. 

4.2 Collider searches for Z' 

One of the main constraint on this model comes from the direct collider searches for a 
Z' boson in the two-lepton decay channel. The best limits have been obtained at the 
LHC by the ATLAS [16] and CMS [17] collaborations for pp collisions at a center-of-mass 
energy of 8 TeV. In [17] limits were obtained with an integrated luminosity of 19.7 fb“^ 
(20.6 fb“^) in the dielectron (dimuon) channel and lead to Mz 2 > 2.57 TeV for Oeq = 0^, 
assuming only SM decay modes. Such limits however depend on the couplings of the Z 2 , 
hence on 9Eg- To reinterpret this limit for any value oi 9Eg, we first simulate Monte Carlo 
signals for Z' production using the same Monte Carlo generator and PDF set as in [16], 
respectively Pythia 8.165 [61] and MSTW2008LO [62], for a large set of 9Eg values. We 
get results compatible with with the ones derived in [63] as well as the one obtained by 
the CMS collaboration [17]. Then we interpolate our limits for any possible choice of 9Eg- 
Note that the coupling of Z 2 to the standard model fermions also weakly depends on az- 
We have checked that this dependence does not modify significantly the Z 2 limits and are 
well below PDF uncertainties [16]. Furthermore in the UMSSM the Z 2 can decay into 
supersymmetric particles, RH neutrinos and Higgs bosons, thus reducing the branching 
ratio into leptons. The limits on the Z 2 mass are therefore weakened [64-67]. To take this 
effect into account we determine in a second step the modified leptonic branching ratio 
for each point in our scan, and re-derive the corresponding limit. 

For any value of 9Eg we restrict the scan to \az\ < 10“^ [49]. In addition, the mixing 
between and Z' can be constrained by the Ap parameter [68]. This observable, which 
measures the deviation of the p-parameter of the standard model from unity, receives 
a specific UMSSM tree-level contribution because Zi is no longer purely the Z^ boson. 
In the limit where M|, ^ which is the case for the TeV scale Z', this new 

contribution reads [68] 

9 ^7 . . 

We compute Ap for each point in the parameter space using a micrOMEGAs routine which 
also contains leading one-loop third generation sfermions and leading two-loop QCD con¬ 
tributions. We impose the upper bound Ap < 8.8 X 10-^ [69]. 

4.3 Collider searches for SUSY particles 

First we impose generic constraints from LEP on neutralinos, charginos, sleptons and 
squarks. For the latter we ignore the possibility of very compressed spectra and use the 
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generic limit at 103 GeV. Lighter compressed squarks can in any case be constrained from 
LHC monophoton searches [70,71] and monojet analyses [72], 

The most powerful and comprehensive constraints on supersymmetric partners have 
been obtained by ATLAS and CMS using the data collected at 7 and 8 TeV. Searches 
were performed for a wide variety of channels and results were presented both in the 
framework of specific models, such as the MSSM, and in the context of simplified model 
spectra (SMS). Here we use the SMS results to find the main constraints on the UMSSM. 
We base our analysis on SModelS vl.0.1 [39,40], a tool designed to decompose the signal 
of an arbitrary BSM model into simplified model topologies and to test it against LHC 
bounds. The version used includes more than 60 SMS results from both ATLAS and 
CMS. 

The input to SModelS are SLHA files that contain the full mass spectrum, decay 
tables as well as production cross sections. The input files, including tree-level production 
cross sections, are generated using micr0MEGAs_4.1.5 [73], for strongly produced particles, 
SModelS then calls NLL-fast [74-80] to compute the k-factor at NLO-I-NLL order. 
Subsequently the code decomposes the full model into simplified model components, and 
calculates the weight (production cross section times branching ratio, a x B) for each 
topology. To limit computing time, topologies with small weights are not considered in 
the decomposition. As a minimum weight we have used a cutoff acut of 0.03 fb. The 
resulting list is then confronted with the SModelS database, for any matching result 
a X B is compared against the experimental upper limit. If no matching results are found 
the point is labeled as not tested. This may happen for several reasons, either all cross 
sections are below acut in which case the decomposition will not return any entries, there 
is no matching simplified model result in the database, or the mass vector of the new 
particles lies outside the experimental grids for all applicable SMS results. Since soft 
decay products cannot be detected, in this analysis we disregard vertices where the mass 
splitting between the mother and daughter SUSY particles is less than 5 GeV. A fully 
invisible decay at the end of a decay chain is compressed, the corresponding mother SUSY 
particle is then treated as an effective LSP. 

Note that topologies that contain long-lived charged particles corresponding to cr > 
10 mm are not tested against SMS results within SModelS. However searches for long- 
lived particles leaving charged tracks in the detector have been performed at the Teva- 
tron [81] and the LHC [82,83] and were interpreted in the context of long-lived charginos 
or in the context of the pMSSM [84]. When the neutralino LSP is dominantly wino, 
typically, the NLSP chargino will be stable at the collider scale. We have therefore con¬ 
sidered the DO and ATLAS upper limits for points with charginos in the mass range 
100 — 300 GeV and 450 — 800 GeV, and decay lengths cr > 10 m and 21 m respectively. 
We have not included the limits from CMS as these cannot be simply reinterpreted for 
direct production of chargino pairs [83]. Long-lived gluinos or squarks are also possible, 
we have not considered these cases since the interpretation of a given experimental anal¬ 
ysis relies on the modeling of R-hadrons, thus introducing large uncertainties. Moreover 
we have not implemented current limits on long-lived staus as these rarely occur in the 
parameter space considered. 
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4.4 Flavour physics 

Indirect constraints coming from the flavour sector, especially those involving 5-Mesons, 
play an important role in defining the allowed parameter space of supersymmetric models, 
e.g. [85-88]. The constraints imposed on the model are listed in Table 2, though we do not 
in general require agreement with the measured value of Aa^. We do however highlight the 
specific regions consistent with the measured value of the muon anomalous magnetic mo¬ 
ment. These mostly correspond to regions with a light LH smuon/sneutrino as mentioned 
in section 2.2. To compute these observables, we have adapted the NMSSMTools routine to 
the UMSSM, for more details see [89]. The most powerful constraints are AM^ and AM^^ 
for small values of tan/3 while ^(5° Wh) and ^(5° ^ are also important 

to constrain some large values of tan/3. We also compute ^(5° but this 

channel does not give additional constraints. Uncertainties coming from CKM matrix 
elements, rare decays, hadronic parameters and theory are taken into account when com¬ 
puting the observables listed in Table 2, see [89]. The most important uncertainties in 
our computation of flavour observables are theoretical (10%) and from the CKM element 
|Ut,| = (42.9 ±2.6) X 10-^ [90]. 


Constraint 

Range 

T^Ur) 

[0.70, 1.58] X10-^ [91] 

^ A37) 

[2.99, 3.87] X10-^ [92] 

^(5° ^ /i+/i-) 

[1.6, 4.2] X10-*^ [93] 

am. 

[17.805, 17.717] ps-i [94] 

am. 

[0.504, 0.516] ps-i [95] 

A (2^ 

[7.73, 42.14] X10-1° [51,52,96] 


Table 2: Flavour constraints used and their allowed ranges which correspond to the 
experimental results (or to the difference between the experimental value and the standard 
model expectation for Ao^) ± 2cr. 


4.5 Dark matter 

The value of the dark matter relic density has recently been measured precisely by the 
Planck collaboration and a combination of Planck power spectra, Planck leasing and 
other external data leads to [97] 


= 0.1188 ±0.0010. (4.2) 

We will impose only the 2cr upper bound from eq. (4.2) on the value of the relic density. 
That is we assume that either there is another component of dark matter or that there 
exists some regeneration mechanism that can bring the dark matter within the range 
favoured by Planck [98,99]. 
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This measurement puts a strong constraint on the parameter space of the UMSSM 
whether the dark matter candidate is the lightest neutralino or the supersymmetric part¬ 
ner of the right-handed neutrino. Since the three RH sneutrinos have the same coupling 
to all other particles in the model we assume for simplicity that the third generation sneu- 
trino is the lightest. In previous studies it was shown that the favoured mass for the RH 
sneutrino LSP was near M 22 / 2 , although much lighter sneutrinos could also be found, es¬ 
pecially near mh^/2 or when coannihilation was present [34]. As in the MSSM the lightest 
neutralino covers a large range of mass, the main new features being the possibility of a 
singlino LSP [100-102] and the possibility for this singlino to have a non-negligeable bino’ 
component. Typical MSSM features can also be observed as the example of wino LSP 
annihilating efficiently into IT’s and strongly degenerate in mass with chargino NLSP. 
However sometimes the mass degeneracy between the NLSP and the LSP can be suffi¬ 
ciently small to give an absolutely stable charged NLSP. When focusing on relic density 
constraints we will systematically discard these configurations. 

One of the strongest constraint on DM arises from direct detection. We implement 
the upper limit from the LUX collaboration [103] taking micrOMEGAs default values for 
the quark coefficients in the nucleons. This upper limit strongly constrains the scenarios 
where the LSP is 0(100 GeV). Another relevant constraint is the one from FermiLAT 
searches for DM annihilation from the dwarf spheroidal satellite galaxies of the Milky 
Way where limits obtained for DM annihilation into bb and can constrain scenarios 

with DM masses below 100 GeV [104]. 


5 Results 


Parameter 

Range 

Parameter 

Range 

^i>rR 

[0, 2] TeV 

p,Mi 

[-2, 2] TeV 

Mz2 

[2.2, 7] TeV 

M2, Ax, At, Ab, Ai 

[-4, 4] TeV 

M[ 

[-20, 20] TeV 

M3 

[0.4, 12] TeV 


[-7r/2, 7r/2] rad 

mp,,miy. 

[0, 4] TeV 

az 

[-10-^ 10-3] rad 

mt 

173.34 ± 1 GeV [105] 


Table 3: Range of the free parameters where concerning the soft mass terms we define 
F G {Q, u, d, L, e}, i G {1,2,3} and j G {1, 2} and where = mt,^. 


After imposing universality for the sfermion masses of the first and second genera¬ 
tion and fixing the trilinear coupling of the first two generation sfermions to 0 GeV, the 
UMSSM features 24 free parameters. The range used for these parameters in the scans 
are listed in Table 3. In addition we have allowed the top mass to vary. We perform a ran¬ 
dom scan over the free parameters and impose first the set of basic constraints described 
from section 4.1 to section 4.3 : the Higgs mass and couplings allowed by HiggsBounds, 
HiggsSignals and our modified NMSSMTools routines, perturbative Yukawas for top and 
bottom quarks, agreement with LEP limits on sparticles and LHG limits on the Z' and 
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finally a neutral LSP. We then include the constraints from 5-pliysics mentioned in sec¬ 
tion 4.4. Another scan is done to highlight the regions of parameter space which give 
sufficient New Physics contribution to Aa^. For this we restrict the soft masses of the 
second generation of sleptons to [0, 2] TeV and we impose all constraints given in sec¬ 
tion 4.4. 



85 90 

m‘;“(GeV) 


0.25 0.3 0.35 0.4 


Figure 3: (a) as a function of the tree-level component of and (b) tan/5 as a 
function of A. For both plots Mz 2 is taken as colour code. 


For all points that satisfy these sets of constraints in both scans, around 4 x 10^, we 
found that the maximum tree-level mass for the Higgs reached only ~ 107 GeV and 
was above the Zi mass only for mixing angles > 2 x 10“^, see figure 3a. Thus a 
contribution from the radiative corrections in the stop/top sector is still required to reach 
a Higgs mass of 125 GeV. Nevertheless the full range of values of tan (5 is allowed. Small 
values of tan /5 > 1 require a large value of A to compensate the small MSSM-like tree-level 
contribution to the light Higgs mass, see figure 3b. This also means that Vg, hence Mz 2 i 
cannot be too large given the range assumed for the /r parameter, see eq. (2.3). Radiative 
corrections from the top/stop sector are expected to be large for tan/5 < 1 since the top 
Yukawa coupling increases as 1/sin/5, which explains why a larger range for A is allowed 
when tan /5 < 1. 

It is well known that large one-loop corrections from the stop sector require heavy 
stops and/or large mixing [8]. The mixing parameter Xt = At — /i/ tan /5 is indeed found 
to be large when Ms = < 1 TeV while heavy stops (associated with large Ms) 

allow no mixing, see figure 4a. The heavier the Z 2 the larger the minimal value for the 
scale Ms where zero mixing is allowed. 

The spectrum for supersymmetric particles differs significantly from the case of the 
MSSM and NMSSM, depending on the choice of V(l)' charges. The lightest stop mass 
can be as light as 300 GeV for Oeq ^ 0.66 (figure 4b), this value corresponds to the 
largest negative contribution to the stop mass from the D-term, see section 2.2. When 
9eq < 0 the lightest stop is at least 670 GeV. Similar values are found for both LH and 
RH up-type squarks, modulo mixing effects. Such light squark masses are well within the 
range of exclusion of LHG searches within the MSSM, hence the need to reinvestigate the 
impact of these searches within the UMSSM discussed in the next section. The dn mass 
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Figure 4: (a) Xt as a function of Ms and as a function of (b) (c) and (d) 

For all plots Mz 2 is taken as colour code. 


receives a large negative F>-term contribution for 9eq = — tan“^(3A/3/5) ~ —1.16. For 
this value it can be as light as allowed by LEP (103 GeV), see figure 4d. For Oeq > 0 the 
RH d-squark is above the TeV scale while the LH one can be light since = rriu^. This 
implies also that a light sbottom, say below 500 GeV, can be found for either value of 
Oeb, see figure 4c. In one case it is mostly LH and in the other RH. Note that an increase 
in the lower limit on the Z' mass will lead to larger squark masses except for the specific 
values of Oeq where one gets a very large T>-term contribution. Finally, the gluino mass 
is determined by M 3 , hence can also be well below the TeV scale. 

The impact of the flavour constraints is best displayed in the tan P — Oeq plane, see 
figure 5. As expected AM^ and AM^ are the most important constraints in our scans and 
exclude a large part of the parameter space when tan/5 < 1, through the charged Higgs 
contribution [89]. The contribution from Double Penguin diagrams to these observables 
enable exclusion of a few scenarios at large tan/5. —)■ and ^sl) are 

important for scenarios at very large tan /5 but they mostly fail to exclude points, especially 
for cases where the mass of heavy neutral and charged MSSM-like Higgs bosons is above 
several TeVs. Finally the New Physics contribution to the deviation of the p-parameter 
from unity exclude only few points, mostly from the sfermion contributions. Actually the 
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pure UMSSM contribution shown in eq. (4.1) can barely reach 10“^ for the allowed values 
for az and Mz 2 and is then negligible. Note that, as we will see in the next section, 
specific regions of the parameter space give large enough contributions to the anomalous 
magnetic moment of the muon. 


10 ^ 


10 


1 


- 1.5 -1 - 0.5 0 0.5 1 1.5 

Be, (fat') 

Figure 5: Points of the scan in the tan/5 - Oeq plane where the colour code shows the 
flavour process that provides the main exclusion. The region that is compatible with Aa^ 
is also displayed. The flavour observables are computed with the NMSSMTools routine 
adapted to the UMSSM. 



5.1 Aa^ 

Special conditions are required to get agreement with the value of Aa^. Indeed the 
discrepancy between the theoretical and experimental value requires a large contribution 
from New Physics. In the UMSSM this comes in particular from diagrams involving 
smuon (LH sneutrino) and neutralino (chargino) exchange. A large UMSSM contribution 
requires either a light smuon/LH sneutrino or an enhanced Yukawa for the muon. The 
latter is found at very large values of tan /5, see figure 6a. A light LH smuon mass arises for 
Oeq = — tan“^(3i/3/5) ~ —1.16 corresponding to a large negative D-term contribution 
as explained in section 2.2. Future collider limits on the Z' mass, say above 5 TeV, will 
severely constrain scenarios for positive values of Oe^ that are in agreement with the latest 
value of Aa^, see figure fib. Note that the distribution of points in the Oe^ — Tnjx^ plane is 
similar to the one found in the general scan where consistency with the muon anomalous 
magnetic moment is not required, except that heavier sleptons are allowed in that case. 

5.2 Dark matter relic abundance 

In this model the LSP can either be a neutralino or a RH sneutrino. The annihilation 
properties of the neutralino LSP are determined by its composition (figure 7). As in the 
NMSSM, the pure bino or singlino LSP is typically overabundant unless it can benefit 
from a resonance enhancement. Note that in this model the Higgs singlet is very heavy 
so that resonant annihilation of a singlino through the Higgs singlet works only for heavy 
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Figure 6: Points allowed by Aa^ in the Oeq — plane, the colour code corresponds to 
(a) different values of tan /3 and (b) • 


singlinos^. The dominantly singlino LSP is found only for masses above 250 GeV. Some 
admixture of a higgsino/wino component or coannihilation processes can however reduce 
the relic density to flh^ 0.1 for any mass. Coannihilation can occur with gluinos or 
other gauginos as well as with sfermions. As in the MSSM the dominantly higgsino or 
wino LSP annihilates very efficienty into gauge boson pairs and therefore leads to an 
under-abundance of dark matter unless the higgsino (wino) LSP mass is roughly above 1 
(1.5) TeV. Note that the B' component of the LSP is never dominant, because the vev 
of the singlet, which mostly drives the mass of the S and the B', eq. (2.15), is always 
above 6 TeV. For \M[\ <C |ns|, S and B' are both shifted towards large masses whereas 
for \M[\ ^ |us| the singlino benefit from a seesaw-type mechanism which allows a singlino 
LSP down to 250 GeV. This close relation between B' and S is illustrated in figure 8. 

We note that the fraction of points that satisfy the 2cr Planck upper bound is much 
higher in the scan where we impose the constraint on Ao^u than in the general scan. The 
main reason is that it is easier to satisfy the relic density upper bound with a bino LSP 
when the sleptons are light. 

Sneutrino dark matter is typically overabundant as sneutrino annihilation channels 
are not very efficient. Agreement with the upper bound set by Planck requires either 
~ mh^/2 or as found in [34]. The latter case requires above the TeV 

scale when considering current limits on the Z' mass, here we consider DM below 2 TeV. 
Annihilation into W or Z\ pairs through Higgs boson exchange was also found to be 
efficient enough for ~ 100 GeV [34]. However this process, which depends mostly on 
the singlet nature of the Higgs boson exchanged, will not give a large enough contribution 
if the lower limit on Mz^ increases as shown in figure 9. Sneutrino LSP masses in the 
range 100 — 1000 GeV are also allowed if some coannihilation mechanism, involving c.g. 
the lightest neutralino or other sfermions, helps reduce the relic abundance. The low 
density of points in this region (see figure 9) reflects the fact that the importance of such 
coannihilation processes require the adjustment of uncorrelated parameters in the model. 

^For an analysis of a scenario with a light singlino DM see [18]. 
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Figure 7: Relic density of B (green), W (red), H (blue) and S (orange) LSP. The 2cr 
upper bound from Planck is shown in grey. 



m^(GeV) 

Figure 8: B' component in the neutralino LSP as a function of its mass with the S 
component in the neutralino LSP as colour code. 


6 Impact of LHC searches for SUSY particles 

After having imposed the basic constraints, flavour constraints and an upper bound on 
the relic density Vth? < 0.1208 (corresponding to the 2cr upper limit of eq. (4.2)), we 
next consider the impact of LHC searches for SUSY particles based on SMS results and 
using SModelS. To analyse the impact of the SMS results we group the points into 
four categories. Points excluded by SModelS are labeled as excluded, points where the 
SMS results apply but the cross section is below the experimental upper limit are labeled 
as not excluded. Points where no SMS result applies, as explained in section 4.3, are 
labeled as not tested. Finally points with long-lived particles cannot yet be tested in 
SModelS. Points that are not excluded are then examined in more details to determine 
the signatures that could best be used to further probe them with upcoming data. We 
divide the study in three steps. First, we consider scenarios with a neutralino LSP and 
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Figure 9: Relic density for Dt-r LSP with Mz^ as colour code. The 2cr upper bound from 
Planck is shown in grey. 


hnd that the most stringent constraints on supersymmetric particles are obtained for 
light gluino or light squarks [6,106,107]. Second, we concentrate on points compatible 
with the measurements of the muon anomalous magnetic moment and that still have a 
neutralino LSP. This dedicated scan provides a significant number of points with light 
sfermions and allows us to ascertain the impact of slepton searches. Finally we investigate 
scenarios with a RH sneutrino LSP, among these we do not characterize the ones that 
are compatible with the muon [g — 2) because of the small number of points involved. 
The possibilities to probe all points with long-lived charginos are considered separately 
in section 6.4 regardless of the dark matter candidate. Our results for the constraints on 
the SUSY spectra are presented in section 7 where we combine all sets. 

6.1 Neutralino LSP 

In most points with a neutralino LSP, the LSP is actually either dominantly wino or 
higgsino, see figure 7. Points with a wino LSP are however mostly not considered in the 
SModelS vl.0.1 analysis because they lead to long-lived charginos. Therefore the most 
common configuration for the supersymmetric spectra relevant for SMS results is one with 
three dominantly higgsino particles with similar masses : the LSP, the second neutralino 
and the lightest chargino. Moreover since the jets/leptons produced in the decay of the 
chargino (second neutralino) to the LSP are too soft to be detected the chargino (second 
neutralino) will often lead to a missing Et (MET) signature. We will see that this has 
important consequences when using the SMS results. In particular hardly any points can 
be excluded from electroweakinos searches as only few can exploit the decay channel into 
real gauge/Higgs boson. Furthermore we do not find constraints from decays into leptons 
via sleptons since sleptons are rarely fight. 

6.1.1 Gluino constraints 

In figure 10 we show points with a neutralino LSP in the LSP and gluino mass plane for 
gluino masses up to 1200 GeV. On the left we show excluded points in red and allowed 
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points in blue, moreover we indicate points with long-lived sparticles that cannot be tested 
in SModelS vl.0.1 in green and points not tested for the other reasons mentioned before 
in grey. The right panel indicates the topology giving the strongest constraint for each 
excluded point. 
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Figure 10: Exclusion with SModelS vl.0.1 in the LSP - gluino mass plane, (a) shows 
whether a point can be tested, and excluded, as well as points which cannot be tested 
because of long-lived sparticles or other reasons, (b) shows the most constraining topology 
for all excluded points. For the most frequently found topologies we specify the associated 
experimental searches : J = [108], f = [109], © = [110], § = [6], 0 = [111], * = [106], o 
= [112], f = [107] and 0 = [113]. 


We find that gluino topologies, basically from gluino decaying into a pair of quarks 
and the LSP through virtual squark exchange, can exclude gluino masses up to 1100 GeV 
[6,107,109]. The exclusions differ from those of a simplified model, since in general there 
are many possible decay channels. The decay branching ratios of the gluino depend 
strongly on the nature of the LSP. For a higgsino LSP, the decay of the gluino via virtual 
stop is dominant because of the stronger coupling which depends on the top mass, the final 
state is thy® (when there is enough phase space) and/or tbx~ where the chargino is treated 
as an effective LSP. The strongest constraints are found when phase space allows only 
the decay into the chargino final state, as there is one dominant decay channel. In other 
scenarios (non-higgsino LSP) there is no such strong preference for one decay channel, and 
the signal cross section will be split up on several simplified model topologies. Moreover 
mixed decays, where each gluino decays into different quark pairs and the LSP occur 
frequently and are not constrained by SMS. Hence the exclusion will be considerably 
weaker than for the pure simplified model exclusion. For many configurations gluinos 
can decay to heavier gauginos yielding topologies with long cascades not yet included in 
SModelS. Moreover each different topology resulting from such processes is typically 
suppressed because of multiple branching fractions. Similarly points with gluino decaying 
via an on-shell sbottoms are not yet included in SModelS while those decaying via an 
on-shell stops can be tested by SMS. However we found that the cross sections are too 
small by two orders of magnitude for these points to be excluded. 
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Points with very light gluinos (below 500 GeV) may remain allowed even for light 
LSP (less than 200 GeV) if the branching ratio g tbx~ is dominant. This is because 
constraints in the region where rrig < 500 GeV available from ATLAS [114] (where the 
chargino is considered degenerate with the LSP) are very weak. This search was also 
considered in CMS [115] but the results are not incorporated in SModelS vl.0.1 as 
digitized data are not yet available. Furthermore results for this topology when the 
chargino is not degenerate with the LSP are only available for one specific mass ratio and 
therefore cannot be used. 

We also found that most points with a light gluino and a dominantly singlino LSP 
feature a very compressed spectrum. This follows from the relic density constraint that 
favours coannihilation as mentioned in section 5.2, thus these points will be hard to 
constrain from SUSY searches for gluinos. 

6.1.2 Squark constraints 

The model can naturally give light squarks, as was shown in section 5. However we 
observe that these are poorly constrained by the SMS limits. We show the excluded vs. 
allowed points as well as the most constraining topology for each point, here in the plane 
of the LSP and the lightest squark mass (including stop and sbottom), see figure 11. 
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Figure 11: Exclusion with SModelS vl.0.1 in the LSP - squark mass plane below 1.5 
TeV, where we select the mass of the lightest squark (including stop, sbottom). (a) has 
same colour code as in figure 10a. (b) shows for all excluded points the most constraining 
topology. For the most frequently found topologies we show the associated experimental 
searches : § = [6], o = [112], 1 = [107], ® = [116], o = [117] and V = [118]. 


A first observation is that U* and 2“*^ generation squark topologies can exclude points 
up to rather high squark masses (about 1200 GeV) in excess of the simplified models 
exclusions. This is expected since a light gluino will enhance the squark production cross 
sections. Note that it was verified in [119] that SMS results can still be safely applied in 
this case. Points along the kinematic edge can in general only be excluded by one heavier 
squark in the point, as such a compressed spectrum cannot be tested by the SMS results. 
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We find however that many points with light squarks remain unconstrained, even 
for large mass differences to the LSP. The first reason for this is that the simplified 
model exclusions depend critically on the assumption that the 8 squarks of the first and 
second generation are nearly degenerate. This is not the case in the UMSSM, where 
because of the new D-term contributions the mass of the RH d-type squarks can differ 
significantly from the other squark masses. Often their masses are not close enough to 
combine the production cross sections before comparing against an upper limit result [72]^. 
We therefore find much weaker exclusions. The second reason is again tied to the nature 
of the LSP. Recall that most points, and in particular the unexcluded ones, feature a 
higgsino LSP, as shown in figure 12, and that important signatures of light squarks with 
a higgsino LSP are not covered by existing SMS results. 
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Figure 12: Higgsino component of the LSP in the LSP - lightest squark mass plane below 
1.5 TeV, for points that cannot be excluded by SModelS vl.0.1. 

To identify the main signatures for squarks that are not covered by SMS results, we 
discuss next the dominant missing topologies, separately for first/second generation and 
third generation squarks. A simplified model topology for which no matching experimental 
interpretation exists is labeled as “missing topology”. 

The notation used for missing topologies keeps track of the branch and vertex struc¬ 
ture. One branch is contained inside brackets, vertices are separated by a comma. Only 
outgoing R-even particles in a given vertex are specified, light quarks and gluons appear as 
jets (denoted by “j”) while third generation quarks are denoted by their name. MET from 
an outgoing DM candidate is always implied, and if no visible R-even particles appear 
in a branch it is denoted as “(inv)”. An example is stop pair production, with ii txi 
in one branch and ij feyf, xt in the other (i,j G {1,2}), illustrated in 

^If several particles (such as squarks of different masses) contribute to the same topology, they will 
be combined if the corresponding masses are found to be compatible. The criterion is based on the 
difference in the experimental upper limits: if they differ less than 20% (while the mass values may differ 
up to 100%), the contributions are merged. For a more detailed explanation see [39]. This is evaluated 
for each experimental result and may hence differ for different experimental analyses considering the 
same topology. Note that despite differences in the upper limits, the contributions of different mass 
configurations may still contribute to the same signal region. Using the appropriate efficiencies for each 
contribution might therefore improve the limits. 
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Figure 13: Sample diagrams illustrating missing topologies and their notation in the 
case of stop pair production followed by an asymmetric decay to a neutralino LSP. Pair 
production, determining the branch structure is shown in red, R-even final state particles 
are indicated in blue. 


figure 13. This topology is denoted as “(t)(6, bF)” if the W is on-shell (figure 13a). In 
scenarios with an off-shell W (as shown in figure 13b) only the decay products will be 
listed, e.g. “{t)(b,jj)” for hadronic W decay (case 1). Finally if the mass gap between 
the chargino and the neutralino is smaller than the limit chosen for mass compression, 

the chargino decay is considered invisible, the topology is then listed as “(t)(6)” (case 2). 
6 

In figure 14a we show the missing topology with the highest cross section for points 
labeled as not tested or not excluded and with light first/second generation squarks. Here 
we select only points where the higgsino fraction in the LSP is greater than 80%. More¬ 
over, to concentrate on topologies derived from squark production, we have removed any 
topology where one branch is fully invisible, thus getting rid of direct higgsino production. 
Indeed in chargino-neutralino production, the neutralino LSP leads to an invisible branch, 
moreover a chargino can also lead to an invisible branch when it is nearly degenerate with 
the neutralino since the soft jets that result from its decay cannot be detected. We further 
remove points in which the dominant missing topology has a weight smaller than 1 fb. 
We find that a main missing topology consists of 4 jets -|- MET deriving from one squark 
decaying to qxi and the other to qx 2 with the neutralino further decaying to the LSP via 
off-shell Zi, giving the additional jets (mostly soft jets). A re-interpretation of the multijet 
analyses for this simplified model could be useful in constraining these points. Note that 
this topology is common in the compressed region where the squark - LSP mass difference 
is small. We also find 3 jets -|- MET topologies, stemming from squark-gluino production 
as described above. These are found mainly when both gluinos and squarks are light 
and the gluinos decay into a squark and a quark. Results for squark-gluino production 

®This notation directly translates to the SModelS nested bracket notation, where nested square 
brackets indicate the branches and vertices. The given example “(t)(6, IF)” is then written as 

[[[t]],[[b],[W]]]. 
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Figure 14: Missing topologies with highest cross section, for higgsino LSP, shown in (a) 
in the m^o — rriq^ plane where gi is the lightest squark of the first and second generations. 
In (b) the mass of the lightest third generation squark is used. For both plots only masses 
below 800 GeV are displayed. Here we have not considered direct higgsino production 
which cannot be tested by 8 TeV LHC results. 


within SMS exist only for almost mass degenerate gluinos and squarks. Note that for 
such points gluino pair production remains unconstrained as the gluino preferably decays 
to on-shell squarks, for which there are no SMS results. Similarly in scenarios where the 
gluino is lighter than the squark, squark pair production remains unconstrained as they 
decay dominantly via gluinos. 

In case of larger squark - LSP mass splittings, we often find gauginos with a mass 
between those of the squark and the higgsinos. In this configuration, the squark can decay 
either to the LSP or to a heavier gaugino, that then decays into the LSP and a gauge boson 
or a Higgs (real or virtual). In particular we find that an important missing topology is 
the one where each pair-produced squarks decays to a different channel, IF)”, but 

“(i; IF)(j) IF)” is dominant in a few points. 

The limits on the third generation squarks are also much weaker than in the simplified 
model. The reason is similar to the one invoked for gluino limits : with the higgsino 
LSP, a stop may decay either to ty® or to bx'^. Therefore, only a fraction of the total 
cross section can be constrained by the simplified model upper limit. Furthermore the 
“mixed” decays, where one of the pair produced stop decays via top and the other via 
bottom cannot be constrained as there are currently no SMS result for this channel. This 
shows up as an important missing topology, “(6)(t)”, in figure 14b. The situation will 
be improved when efficiency maps will be incorporated into SModelS [120]. When the 
mass splitting between the stop and the LSP is below the top mass, the main missing 
topology is rather associated with sbottom pair production with one sbottom decaying to 
bxi and the other to bx 2 , followed by X 2 QQXi via an off shell Zi leading to ‘^{b){b, jjy\ 
Similarly “(t)(t,ii)” appears at large mass splittings. An important missing topology is 
the one associated with chargino pair production with charginos decaying to the LSP and 
jets or leptons via a virtual IF, “(ii)(i^e)”. We further find a few points where direct 
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production of heavy charginos, decaying via W, gives the dominant missing topology 

‘Xw)(wy’. 

Note that listing missing topologies with the largest cross section can sometimes be 
misleading as the background was not taken into consideration. It is certainly possible 
that a signature with a smaller cross section gives a better signal to background ratio. 
Examples are leptonic vs. hadronic decays of the W, or decays into b-quark as compared 
to decay into light jets. 

6.2 Neutralino LSP : Aa^ and slepton constraints 
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Figure 15: Exclusion with SModelS vl.0.1 in the m^o — rrij^^ plane below 500 GeV. (a) 
has same colour code as in figure 10a. (b) is showing for all excluded points the most 
constraining topology. For the most frequently found topology we show the associated 
experimental searches : 0 = [121] and 0 = [5]. 


We have separately studied points where the muon anomalous magnetic moment con¬ 
straint is fulfilled. Because of the light smuons (and selectrons) additional LHC constraints 
from slepton SMS topologies become relevant. These constraints played a marginal role 
in the general scan considering the small fraction of points with light sleptons. We show 
the exclusions in figure 15. Excluded points are found mainly in a small region in the 
mass plane, for light smuon masses between 250 and 380 GeV. These exclusions 
are slightly weaker than the ones obtained in ATLAS and CMS [5, 121] which assume 
all sleptons decay into while here sleptons can also decay into . Moreover for 
weakly interacting particles we only compute the production cross section at leading-order 
while SMS include NLO cross sections. A single point is excluded by the slepton SMS 
result despite a very small mass difference between smuon and LSP. However, in this case 
it is not actually the slepton production that is being constrained, but pair produced 
charginos, each of them decaying to a left handed sneutrino which then decays invisibly 
to the neutralino. The signature is identical to that of slepton pair production, giving a 2 
lepton and missing energy final state and was discussed in [122]. Other exclusions come 
from gluino and squark topologies, as described in the previous section. 
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6.3 RH sneutrino LSP 


In the case of a sneutrino LSP we have to bear in mind that since these sneutrinos are RH 
all decays of heavier sparticles must proceed via a neutralino. When the neutralino is the 
NLSP it decays invisibly into or prUr, therefore signatures resemble those associated 
with a neutralino LSP. When decays through an on-shell neutralino are not allowed, we 
effectively find additional neutrinos in the decay vertex to sneutrino, for example in the 
squark decay q —?• qurUr. The signature is essentially the same as for a squark decaying 
into a neutralino LSP since the neutrino will contribute only to the MET, but the event 
kinematics can be changed due to the additional invisible particle in the vertex. This 
issue remains to be investigated and these signatures are not treated in SModelS vl.0.1. 
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Figure 16: Exclusion with SModelS vl.0.1 for points with RH sneutrino LSP (a) in the 
neutralino - gluino mass plane and (b) in the neutralino - lightest squark mass plane. 
Points with very heavy neutralinos and squarks are not displayed. 


Figure 16 is showing points with a RH sneutrino LSP in the Xi ~ 9 Xi ~ 9i 
mass planes. One striking feature is that in this scenario there are many points with 
long-lived gluinos or squarks. Those are mainly found when the lightest neutralino is 
heavier than the gluino or squark since decays into RH sneutrino LSP can only proceed 
via some virtual sparticle and are hence suppressed. Among the points that can be tested, 
only a small number can actually be excluded. The exclusion channels are similar to the 
ones for the neutralino LSP discussed in sections 6.1.1 and 6.1.2 and involve a decay of 
a gluino or squark through a neutralino which further decays into the LSP. We find no 
exclusion from electroweak production. It is therefore instructive to consider the missing 
topologies. To clarify again the notation used for missing topologies, we show in figure 17 
the case of chargino-neutralino production for a sneutrino LSP. The neutralino decays 
to a neutrino and a sneutrino, making this branch entirely invisible, hence indicated by 
“(inv)”. The signature of the chargino decay will depend on whether the intermediate 
neutralino is on-shell or not, indicated by cases 1 and 2. If the neutralino is on-shell its 
decay will be invisible and it can be considered as an effective LSP, yielding the topology 
“(inv)(IT)”. If on the other hand the decay to on-shell neutralino is not possible, the 
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chargino will effectively decay directly as yf the topology is then described 

as “(inv)(hhz^)”. Recall that in the SModelS nested bracket notation, this topology is 
denoted by [ [] , [ [W, nu] ] ]. 



case 1 : (inv)(tP) 
case 2 : (inv)(tP;^) 

Figure 17: Diagram illustrating missing topologies and their notation in the case of 
chargino-neutralino production where the sneutrino is the LSP. Pair production, deter¬ 
mining the branch structure is shown in red, R-even final state particles are indicated in 
blue and the invisible branch is represented in green. 

We show (for all not excluded or not tested points) the missing topology with the high¬ 
est cross section, selecting only the five most frequent ones in each plane, see figure 18. 
At low neutralino masses topologies from neutralino-chargino production are often dom¬ 
inant, with the charginos decaying either directly to WprPr, “(inv)(fFz 2 )”, or via W*Xi, 
“(inv)(jj)”. In both scenarios the neutralino decay is invisible. Note that for the miss¬ 
ing topologies we do not distinguish between LH or RH neutrinos. The SMS limits on 
chargino-neutralino production with W^*'^ final state cannot be applied for either topol¬ 
ogy. The reason is that SMS results assume that the process involves one of the heavier 
neutralinos which then decays via a gauge boson and the LSP, whereas here only the 
chargino decays into visible particles. Moreover, in the first case, there is an additional 
neutrino in the decay. In the second case the decay products of the W* are very soft 
because of the degeneracy between the lightest chargino and neutralino. Pair produced 
charginos decaying to Wxi also provide an important topology, “(IF)(IF)”. Both the 
lightest and heaviest chargino can contribute to this topology. A similar topology with 
off-shell IF’s also occurs although it is suppressed by the hadronic branching ratio. Note 
that current SMS results for this topology only give weak constraints and are not yet 
included in SModelS. 

Finally topologies associated with squark pair production are also frequently dominant, 
for light squarks and heavier neutralino we find squarks decaying directly to the right 
handed sneutrino, q qprPr with either a light quark or a b-quark, corresponding to the 
topologies “{pj){pjy^ and ^\pb){pby' in figure 18. Missing topologies involving gluinos 
are similar to the ones in the neutralino LSP case, note however that when lighter than 
the gluino is likely to be long-lived, or otherwise to decay via 4-body, g jjPRPR- 
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Figure 18: Missing topologies with the highest cross section in the neutralino - light¬ 
est squark mass plane, for points with a RH sneutrino LSP that are not excluded by 
SModelS vl.0.1 and that do not involve long-lived sparticles. 


6.4 Long-lived charged NLSP 

The DO collaboration has searched for pair produced long-lived charginos [81], putting 
upper limits on the production cross section for chargino masses between 100 and 300 
GeV. Since experimental limits are given separately for wino and higgsino-like chargino, 
we use the relevant result and in case of large mixing {i.e. wino fraction in yf between 
0.3 and 0.7) we apply the more conservative limit. Note that the limit is only marginally 
different in the two cases. Results are shown in figure 19a. We find that long-lived 
charginos lighter than about 230 GeV are excluded. 

In addition, the ATLAS collaboration has searched for long-lived charginos from either 
pair production of charginos or chargino-neutralino production [83] , yielding upper limits 
on the combined cross section for chargino masses between 450 and 800 GeV. We have 
checked that the less constrained chargino-neutralino contribution is never larger than in 
the scenario considered by ATLAS, thus ensuring that the application of the upper limit 
is always conservative. Note that in addition to chargino pair production we generally 
consider only yf y? production, except when this is essentially zero then we include also 
xfx 2 production. This may occur if the LSP is bino or singlino are degenerate in mass 
with the chargino^. Results are shown in figure 19b. We find that even at low masses some 
points cannot be excluded, because interference between light squark exchange diagrams 
lead to small production cross sections. However, a large number of points, with chargino 
masses up to about 650 GeV, can be excluded. Note that in both cases we have used 
linear interpolation between the given data points. We expect that smaller masses (below 
450 GeV) should be excluded as well, but existing searches in that mass range consider 
long-lived staus (ATLAS) or long-lived leptons (neutral under SU{2)l, see [83]) and were 
not applicable here. 

^This degeneracy can follow from imposing the relic density upper limit which in this case will be 
satisfied because of the contribution from efficient coannihilation channels involving the chargino and 
heavier neutralinos [123]. 
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(a) 



Figure 19: Points tested by searches for long-lived charginos: (a) the chargino pair pro¬ 
duction cross section and the corresponding upper limits from DO (b) sum of chargino 
pair and chargino neutralino production cross sections and the corresponding upper limit 
from ATLAS. The colour code indicates the wino fraction in xt- 


Finally we point out the potential of such a search at 13 TeV. In figure 20, the cross 
section for pair production of charginos with decay lengths cr > 10 mm is displayed. Here 
all points that have not yet been excluded are shown. We find that about one order of 
magnitude improvement over the current limit would allow to probe a large fraction of 
the points with a long-lived chargino below the TeV scale. Note that in this figure we 
have included long-lived charginos decaying either inside or outside the detectors. Each 
category includes a significant number of points. Therefore both types of searches could 
be used to test the model further. 


7 Summary after LHC constraints 

7.1 Exclusion potential of current LHC searches on the UMSSM 

To summarize the impact of the LHC constraints on the sfermion spectrum we display in 
figure 21 the excluded/non-excluded points in the plane Oeq — for / e {ii,bi,dR} as 
well ss f = fiL for the sample where the muon anomalous magnetic moment constraint is 
imposed. Among the non excluded points those that satisfy all constraints have a different 
colour code than those that are associated with a long-lived NLSP or that are not tested 
by SModelS vl.0.1. In all cases the excluded points are scattered and represent only 
a fraction of all points. It should be stressed again that many scenarios with squark 
masses well below 1 TeV are allowed. When the agreement with Aa^ is not required 
we found that 45% (41%) of the points that were confronted with the LHC limits had 
a long-lived sparticle in the case of a neutralino (RH sneutrino) LSP, 16% (17%) were 
tested by SModelS of which 10% (11%) were excluded. The remainder of the points 
was not testable by SModelS either because of too low cross sections or lack of SMS 
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Figure 20: Chargino pair production cross section at 13 TeV, for unexcluded points with 
long-lived charginos where the colour code indicates either (a) the wino fraction in xf or 
(b) CT < 100 m. 


result. We additionally found that 42% (24%) of the sample with long-lived NLSP were 
excluded by long-lived chargino searches. In the case where the muon anomalous magnetic 
moment constraint is required and for a neutralino LSP the amount of points tested by 
SModelS is larger (34%, out of which 11% are excluded), whereas the fraction of points 
with long-lived sparticles is smaller (30%, out of which 44% can be excluded). Extending 
these searches for long-lived charginos to the full mass range would therefore clearly 
provide a powerful probe of the model. Moreover, when the RH sneutrino is the LSP a 
large fraction of the points involves long-lived gluinos and squark. These scenarios could 
test the model further, but require reliable limits on R-hadrons. Note that to facilitate 
the interpretation of limits on long-lived charginos, it would be useful if limits on the 
direct chargino pair and neutralino-chargino production were separately provided by the 
experimentalists. 

Many of the points that are in agreement with the measured value of the muon anoma¬ 
lous magnetic moment, even those associated with a very light smuon cannot be com¬ 
pletely excluded, see figure 21d. The LHC13TeV with higher luminosity will allow to 
extend the reach for smuons in the conventional lepton -|- MET channel. Unfortunately 
the light charginos that are present in this case cannot be probed easily as once again 
they are often dominantly higgsino hence almost degenerate with the LSP. 

7.2 Suggestions for future LHC searches 

Conventional searches for a new Z' provide the most distinctive signature of the UMSSM. 
In addition we have pointed out in previous sections many additional signatures that 
are still unconstrained by current SMS results. Here we summarize the main missing 
topologies found for each scenario. 

As expected, the most distinctive SUSY signatures in the UMSSM are found in the 
case of a RH sneutrino LSP. We have already stressed that long-lived gluinos or squarks 
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Figure 21: Points in the 9eq — rnj plane for (a) / = ti, (b) / = 6i, (c) / = and 
(d) / = JIl. For all plots < 1.5 TeV and the colours correspond to not tested 
(grey), long-lived (green), not excluded (blue) and excluded (red) configurations using 
SModelS vl.0.1 and the searches for long-lived charginos. For (d) only points satisfying 
Ao^ are represented. 


are fairly common in such scenarios and could therefore provide further constraints on 
the model. We have also found that the following topologies could be used to probe the 
model either by reinterpreting current LHC data or by exploiting Run II data. 

• mono-bF, “(inv)(lFi^)”, from chargino neutralino production with —>■ vrVr and 

Xi —t W'^vrVr. The single W can be energetic enough to lead to visible decay 
products, leptons or jets as long as there is a large mass difference between the 
chargino and sneutrino LSP. Such a topology occurs also for a neutralino LSP (as 
in the MSSM) but only when there is a large — Xi mass splitting to allow 
xf —>■ W^Xi, which is not the most common configuration after imposing DM 
constraints. 

• dijets -|- MET, “(i^j)(z^j)” or bb + MET, “(z/6)(z/6)”, from squark pair production 
with q qi'Ri'R where q here stands for either light jets or b-jets. This occurs 
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when the squark is lighter than all neutralinos and therefore has to decay directly 
to the sneutrino LSP. Such a configuration is clearly only possible in a model with 
a sneutrino LSP. The dijet + MET signature is of course common to squark pair 
production in the MSSM, however it remains to be seen how the additional vr in 
the decay will affect the efficiencies, hence could lead to different excfusions than in 
the case of the neutralino LSP. 

Other important missing topologies include dijets + MET, “(inv)(ji' j)”, from chargino 
neutralino production with the same decays as the mono-lT above except that the W is 
off-shell leading to soft final states as well as W pairs -|- MET, “(1T)(1T)”, from chargino 
pair production with Xi These topologies also arise in the MSSM 

with a neutralino LSP and are poorly constrained from searches at Run 1 partly due to the 
small production cross section. The situation should however improve after accumulating 
more data in Run 11. 

When the neutralino is the LSP, most SUSY signatures are the same as found in 
the MSSM. However we stress that having imposed only the upper limit on the dark 
matter relic density, most of our scenario have a wino/higgsino-like LSP. Thus, the SUSY 
signatures can differ from the bino LSP assumed in severaf SMS resufts. In particular, the 
chargino decay can be invisible and this will have an impact on many SUSY searches. In 
addition this implies that a significant fraction of the scenarios have a long-lived chargino 
and/or neutralino, hence the importance of searches for stable charged particles at collider 
scale and for displaced vertices. 

Many of the topologies that could not be constrained by current SMS results are 
associated with asymmetric decays, that is the pair produced particles have two different 
decay chains whereas most SMS resufts assume identicaf decays for both particles. We 
emphasize here the missing topologies for the case of the higgsino LSP since it is hard to 
probe. 

• 3 jets -|- MET, ^‘{j){j,jy\ from gluino-squark production with g ^ qq and q —?• qx^. 

Current SMS interpretations exist only for scenarios where the gluino and squarks 
are almost mass degenerate, and both decay directly to jets and LSP. Similarly the 
topology 4 jets -|- MET, j)’\ from gluino pair production with g ^ qq 

and q —>■ gy? arises from a process with a large production cross section that is 
not constrained by SMS results since the gluino decays via on-shell first or second 
generation squarks. Note that both these topologies are of special interest in the 
UMSSM where the limits on light squarks from direct squark production are much 
weaker because the squarks are not necessarily all degenerate. 

• bt + MET, “(6)(t)” , from stop (sbottom) pair production with asymmetric decays, 

i and i bx^ (b —?• 6%° and b tyj") when the chargino is nearly degenerate 

with the LSP. This signature is a generic feature of models with wino/higgsino LSP 
and light third generation squarks [112,124]. 

• 4 jets -|- MET, “(i)(i,jj)”, from squark pair production with asymmetric decays. 
Here one squark decays directly to the LSP, q gy? while the other decays via 
heavier neutralino, g —>■ gy® and x^ ^ Z^Xv A re-interpretation of the multi-jet 
analysis to study the effect of the soft jets from the virtual Zi on the efficiency 
woufd be useful. 
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• 26 + 2 jets + MET, “(6)(6, jj)”, from sbottom pair production with asymmetric 
decays, same as above. Note that this topology is found for small mass difference 
between the sbottom and the LSP. Similarly the 26 + 2 jets, “(6)(6, jj)”, from stop 
pair production is also a missing topology. 

• 2 jets + W + MET, “{j){j,W)”, from squark pair production with asymmetric 
decays. One squark decays to LSP q —> qxi while the other decays qu —>• qdXi with 
Xi —t W~^Xi or g ^ with X 2 ^~xt when the chargino decays invisibly. We 
find this when the mass splitting between the squark and the LSP is large. Note 
that typically there would be similar channels with Zi or hi in the final state instead 
of a IT reducing the cross section for each single channel. 

• 2 jets + ITIT + MET, from squark pair production, here both 

squarks can decay to chargino, followed by xf —?• IT^y? above. This channel 
has been considered at the LHC but is not yet included in the SModelS database, 
moreover results are available only for specific mass relations. It would be preferable 
to provide SMS results that allow for interpolation over wide range of masses in 
different mass planes. 

In addition the signature ITIT + MET from chargino pair production will be useful 
in constraining the model, although it currently gives only weak limits [121]. Similarly 
the signature 2 jets + lepton + MET from charginos decaying via virtual IT’s is often 
found, current data do not put useful constraints but these searches should be improved 
in the next Run. 


8 Couplings and signal strengths for the Higgses 

In the UMSSM a lightest Higgs scalar with a mass of 125 GeV is easily found. Typically 
this lightest scalar is doublet-like and behaves roughly as the SM Higgs. Measurements of 
the Higgs couplings at the LHC Run H could therefore provide additional probes of the 
model. Eor all points of the UMSSM scan that successfully pass all collider constraints 
we have computed the reduced couplings of the 125 GeV Higgs. The reduced couplings 
are defined as scaling factors of the couplings in the UMSSM relative to their SM coun¬ 
terparts. We find that the hiW~^W~{hiZiZi) couplings deviate by at most 1% from the 
SM couplings while there is more room for deviations in the quark couplings. The hibb 
reduced coupling (Cb) can be as large as 1.2 for large values of tan/5 while the hitt cou¬ 
pling {Ct) can be suppressed by at most 5% for low values of tan /5, see figure 22a. Note 
that the couplings are generation universal. Modifications of the quark couplings induce 
a correction to the loop-induced couplings of the Higgs to gluons {Cg) and photons {C-y). 
In particular since the top quark gives the largest contribution to Cg in the SM, we expect 
a reduction in Cg. Eurthermore, this should be correlated with a mild increase in Cj as 
observed in figure 22b. Supersymmetric particles can also contribute to the loop-induced 
coupling, for example light squarks can lead to > 1 although the effect is again below 
5%. Light sleptons and charginos will only contribute to C^. Again the effect typically 
does not exceed 5%. 

The effect on the Higgs signal strength can be much larger than on the reduced cou¬ 
pling. The signal strength in one channel is defined as the production cross section times 
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(a) (b) 

Figure 22: Reduced couplings of the light Higgs for all points satisfying collider and DM 
relic density constraints : (a) Cb vs. Ct with tan (3 as colour code, (b) Cg vs C-y with Ms 
as colour code. 


branching ratio in the UMSSM relative to the SM expectation for a Higgs of the same 
mass. An increase in the total width, through an increase of the dominant bb partial 
width, will lead to a reduced branching ratio, hence to a reduced signal strength, in all 
other decay channels. Furthermore, when new decay modes are possible (here it means 
invisible decays into the LSP) the total width of the Higgs increases, thus reducing the 
signal strengths in all channels. For example the signal strength for the two-photon mode 
in gluon fusion can be reduced by 25% as compared to the SM expectation, see 

figure 23a. Because this large reduction comes from the total width we expect it to be 
completely correlated with the signal strength in the W fusion mode. A comparison with 
the signal strengths for the bb mode, figure 23b, clearly shows that this reduction can be 
correlated with the one in the bb channel (when the invisible width is large) or with an 
increase in the signal strength in the bb channel when Cb > 1- The invisible width of 
hi is found to be below 25%. Recall that current limits from direct searches are 58% in 
CMS [125] while preliminary results from ATLAS in the vector boson fusion mode set the 
limit at 29% [126]. A stronger limit of 12% is obtained from global fits to the Higgs [127], 
however the latter applies only when all Higgs couplings are SM-like. In future runs, it is 
expected that the LHC could probe directly an invisible width of 17% [128]. 

Other probes of the Higgs sector can be performed by searching for the heavy Higgses 
at the LHC. After applying all constraints described in section 4, which in particular 
include heavy Higgs searches at LHC8TeV, we find that the lowest allowed value for the 
mass of ^2 is around 340 GeV and can reach several TeV’s. Below the TeV scale, the 
pseudoscalar is typically nearly degenerate with the doublet-like ^2 and the value of tan (3 
ranges from 2-40 with a large fraction of the points with tan/5 < 10 because of flavour 
and direct search limits. 

To compare with the recently released limits on searches for heavy Higgs in the W~^W~ 
channel we have computed the signal strengths for h 2 W~^W~ in both the VBF and 
gluon fusion mode. We expect this signal strength to be quite low (as we have argued 
above the coupling hiW~^W~ is SM-like). This means that in the decoupling limit the 
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Figure 23: Signal strengths of the light Higgs into (a) 77 and (b) bb in the Vector boson 
fusion (VBF) vs. gluon fusion for all points satisfying collider and DM relic density 
constraints. For both plots is taken as colour code. 


h2W^W~ coupling is suppressed, cos(q; — /5) ~ 0 in the MSSM notation. Indeed we hnd 
that the signal strength is suppressed in the gluon fusion channel, fi^{h 2 ) < 0.03, due to 
the small branching into gauge boson final state and obviously even more so in the VBF 
production mode where the signal strength is well below 10“^, see figure 24a. Thus h 2 
easily escapes current limits. Note furthermore that the largest signal strengths are found 
for low values of tan /3 and for /12 much below the TeV scale, a region where potentially 
the tt channel offers a better probe, as discussed below. 

In the sub-TeV region, preferred decay channels of h 2 are usually in the bb (r+r") 
final states for moderate to large values of tan/5, as in the MSSM. However, for low 
values of tan /5, h 2 can decay exclusively into ti, see figure 24b. Moreover decays into the 
lightest Higgs can also be large (as much as 50% when m /12 < 360 GeV) but drop rapidly 
reaching at most 10% when m /12 > 460 GeV. In the MSSM it was shown that searches 
for heavy Higgs in the tt [hh] channel offer good discovery potential at LHC13TeV for 
small values of tan/5, when rtih^ < 1(0.5) TeV [129]. Hence, such searches should also 
probe of the UMSSM model further. However, decays of h 2 into supersymmetric particles 
can affect the main SM particle signatures. In particular decays into electroweakinos can 
reach 84% ( 86 %) for the neutralino (RH sneutrino) LSP scenarios, while the invisible 
decay of ^2 into the neutralino LSP reaches at most 15%. Large branching fractions 
into electroweakinos are expected when the kinematically accessible states have a large 
higgsino/gaugino component, hence when yt, M 2 are small. These decay modes could 
therefore provide additional search channels for a second Higgs, see e.g. [130]. For a Higgs 
below the TeV scale, the decays into sfermions are generally kinematically forbidden. 
When they are allowed the branchings never reach the percent level and are therefore 
negligible. 
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Figure 24: For all points satisfying collider and DM relic density constraints are displayed 
(a) the signal strengths for /12 in gluon fusion and VBF modes in the W~^W~ final state, 
with 771/12 < 1 TeV as colour code, (b) <^(/i 2 —t it) in the tan/5 — rriAO plane with 
rriAO < 1 TeV. 


9 Dark matter probes 

The correlation between LHC constraints and DM relic abundance on both the neutralino 
and RH sneutrino LSP scenarios is summarized in figure 25. Clearly there is a strong 
preference for a RH sneutrino around 60 GeV, see figure 25a. Moreover for the neutralino 
LSP case, the wino scenario (the lower branch in figure 25b) is strongly constrained by 
searches for long-lived charginos. We now consider the predictions for DM observables. 
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Figure 25: Relic density for (a) i^t-r or (b) Xi LSP with the same colour code as in 
figure 21 except for the configurations excuded by searches for long-lived charginos which 
are highlighted in pink. 


Direct searches for dark matter through their scattering on nuclei in a large detector 
provide a complementary method to probe supersymmetric dark matter. When examining 
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the predictions for dark matter searches we use a rescaling factor to take into account cases 
where the LSP constitutes only a small fraction of the dark matter. The 2g deviation 
from the central value measured by Planck is = 0.1168, hence we define the rescaling 
factor ^ as 


for < 0.1168, 

1 for G [0.1168,0.1208]. 


(9.1) 



Figure 26: Rescaled direct detection cross section for (a) u-j-r or (b) y? LSP with the same 
colour code as in figure 25. LUX exclusion (dark beige), projections from future large 
detectors as well as the neutrino background are also displayed. 


Figure 26a shows that most of the RH sneutrino LSP points with a mass around 
100 GeV are excluded by LUX. Those with a mass near 60 GeV escape the LUX upper 
limit but are generally within the reach of the future XenonlT detector. Other RH 
sneutrinos, which as we have argued before benefit from coannihilation and are therefore 
associated with a compressed spectrum, are safely below current exclusions. In some 
cases the predicted cross section is even below that of the coherent neutrino background 
and can therefore never be probed by direct detection. The scenarios with a neutralino 
LSP are hardly probed by LUX, see figure 26b, only some of the mixed bino/higgsino 
points are excluded. The XenonlT will be able to probe many more points, although a 
large fraction is beyond the reach of even a 10 ton detector, if not below the coherent 
neutrino background. These points are dominantly wino (hence labelled as long-lived) 
or singlino LSP. It is interesting to note that many of the points that are out of reach 
of direct detection detectors are associated with long-lived sparticles. To illustrate the 
complementarity with collider searches, we show in figure 27a the points with a long- 
lived chargino which could be probed at LHC Run H, that is the points in figure 20 for 
which the cross section for chargino pair production is above 0.1 fb. Clearly, many of 
the points with charginos stable at the collider scale have a direct detection cross-section 
below the reach of XenonlT and even below the neutrino background. Note that the 
lowest value for the direct detection is about four orders of magnitude below the neutrino 
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background (not shown in the figure). It should also be emphasized that many points 
with a chargino lifetime that leads to displaced vertices (in blue and green in figure 27a) 
are also beyond the reach of ton-scale detectors, hence we stress again the importance 
of probing these signatures at colliders. A quite different conclusion would be reached if 
we set the rescaling factor ^ = 1, that is assuming some regeneration mechanism for the 
neutralino LSP. As shown in figure 27b, some of the mixed wino points are not allowed 
by LUX and a large fraction are within the reach of XenonlT. However even with these 
optimistic assumptions we find a few scenarios with a cross section below that of the 
coherent neutrino background that will never be tested. 



Figure 27: (a) Rescaled direct detection cross section with a i>rR or a Xi LSP for the same 
set of points as in figure 20 and after applying the LUX bound. Chargino decay lengths 
cr > 20 m (red), between 20 m and 1 m (green) and below 1 m (blue) are represented, 
(b) Direct detection cross section without rescaling for B (green), W (red), H (blue) and 
S (orange) LSP, for all points satisfying collider and DM relic density constraints. 


We also explore the possibility to probe DM with indirect detection, in particular 
using the limits obtained from observations of photons from dwarf spheroidal galaxies in 
the Milky Way by FermiLAT. For this we again rescale the cross section for points where 
the RH sneutrino or neutralino LSP cannot explain all the cosmologically measured dark 
matter. This means introducing a suppression by a factor We find that only a few 
points with a RH sneutrino LSP with a mass near 60 GeV are excluded by the FermiLAT 
limit from the bb channel, see figure 28a. Some of these points were also excluded by 
LUX, however the predicted cross sections for most of the points are suppressed by at 
least two orders of magnitude as compared to current limits. For heavier LSP’s the 
preferred annihilation channel is into W pairs. After applying the rescaling factor no 
exclusion can be obtained. Again, the predicted cross section is generally two orders of 
magnitude below the current limit, see figure 28b. A quite different conclusion is reached 
if one does not apply the rescaling factor, then all winos with a mass below 500 GeV are 
excluded by FermiLAT as well as most of the higgsino LSP’s with a mass below 200 GeV, 
see figure 29. This is not a specific feature of the UMSSM and was already observed in 
the MSSM both for photons [131,132] and also antiprotons [133]. Note that the singlino 
LSP scenarios cannot be probed in this channel even without the rescaling. 
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(a) (b) 


Figure 28: Rescaled annihilation cross section into (a) bb (b) W~^W~, for i>rR (green, red 
when excluded by LUX) or Xi (blue) LSP, for all points satisfying collider and DM relic 
density constraints. In both cases the FermiLAT limits are displayed. 



Figure 29: Annihilation cross section into W~^W~ without rescaling for B (green), W 
(red), H (blue) and S (orange) LSP, for all points satisfying collider and DM relic density 
constraints. 


10 Conclusion 

We have reexamined the viability of the UMSSM model to describe physics beyond the 
standard model and dark matter after the results of the LHC Run I on the Higgs, on 
flavour observables and on new particle searches. We found compatibility with all latest 
experimental results for large regions of the parameter space for both a neutralino or RH 
sneutrino DM and explored potential future probes of the model at Run H of the LHC 
and in direct detection. Imposing only the upper bound on the relic density favors either 
sneutrino DM near 60 GeV or neutralino DM with a large wino and/or higgsino component 
- hence associated with either a long-lived chargino or a chargino which decays primarily 
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into virtual W and the LSP. The latter feature is common also to the MSSM and has 
important consequences both for squark and gluino searches as well as for electroweakino 
searches at the LHC, since the chargino will either have a missing energy signature or 
yield soft decay products. 

The most crucial test of the model in Run II of the LHC is the search for a Z' gauge 
boson, the most prominent characteristic of the model. Another landmark of the model 
is the possibility of a RH sneutrino LSP. The collider signatures often resemble those of 
the MSSM but it remains to be seen whether the additional neutrino that appears in the 
decays affects the kinematics of the process, for example this is the case for the dijet + 
MET signature from squark production when the squark can only decay to the sneutrino 
LSP. We have also found a IT+invisible channel that could lead to a single energetic 
lepton from chargino/neutralino production. Due to the presence of the sneutrino LSP, 
the W can be energetic enough to lead to visible decays even when the chargino and 
lightest neutralino are nearly degenerate, thus this is quite a distinct signature from the 
MSSM. Another striking feature of the sneutrino LSP scenario is the possibility to hnd 
long lived squarks and gluinos. 

The main characteristic of the neutralino LSP scenarios is associated with the nature 
of the LSP, we found usually a dominanty wino or higgsino LSP. This implies that the 
chargino is often nearly degenerate with the LSP hence leads to invisible decay products 
or that the chargino is stable at the collider scale, hence the importance of the searches 
for long-lived/stable particles. This is not a unique feature of the model as long-lived 
charginos can also be found in the MSSM. We also found a number of topologies that 
could not be constrained by current SMS results, including 2 , 3 or 4 jets -|- MET topologies 
with either light jets or third generation quarks. The main reason is that these topologies 
are associated with asymmetric decays, that is each of the pair produced particles has 
a different decay chain whereas most SMS results assume symmetric decays. The most 
striking example is a bt + MET topology that comes from stop pair production where 
one of the stop decays into a quark and the LSP and the other to a quark and an invisible 
chargino [134]. Note however that a recent result by CMS [ 112 ] gives limits on stop pair 
production for this topology independently of the relative branching fractions of each stop. 
It can be expected that this result will exclude most points with the bt + MET missing 
topology in scenarios where the neutralino is lighter than ~ 200 GeV and the stop mass is 
below ~ 700 GeV. Other examples include decays of squarks through a heavier neutralino. 
Furthermore we have stressed that there is a nice complementarity between the searches 
for stable charginos at colliders and the direct searches for DM. A significant number of 
the points which are below the reach of ton-scale detectors have a chargino which is either 
collider stable or lead to displaced vertices. Similarly the singlino LSP, which occurs for 
a small fraction of the points, can hardly be probed by DM direct detection. 

Another feature of the UMSM is the split u-type/d-type RH squark masses which are 
found for specific choices of the G(l)' charge. Despite conventional decays, these squarks 
are harder to detect. They do not beneht from contributions of all flavours of squarks 
and therefore the production cross section is lower. 

As concerns the Higgs sector, the model predicts mostly a SM-like Higgs, although 
deviations up to 25% can be observed in the signal-strengths for either 77 or fermionic 
hnal states. Moreover an invisible branching fraction of the Higgs up to 25% can be found. 
Over some of the parameter space the second Higgs lies below the TeV scale and can be 
probed at LHG13TeV in the standard T+r“ mode relevant at large tan (3 but also in the 
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tt or hh mode at small values of tan/3. It remains to be seen to which extent the SUSY 
decay modes can be exploited. Finally we should stress that the lowest values of tan/3 
lead to quite enhanced rates for AMg. Refining the constraints on the CKM elements 
which is one of the important source of uncertainties implies that this observable would 
strongly constrain the low tan /3 scenarios. 
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Appendices 

A Radiative corrections in the Higgs sector 

To introduce radiative corrections in a gauge invariant manner we use an effective La- 
grangian, 


-^e// = Ai|77,|V2 + A2|77„|V2 + \^\Ha\^\H^\^ 

+ A4|/7„ • ■ Hdf + • HdD/2 

+ iXe\Hd\^ + A7|77„P)((i7„ • H^) + {H^ • Ha)*) (A.l) 

+ Xs\Hd\^\S\^ + X9 \Hu\^\S\^ 

+ Xs{SH^-Hd + S*{Hu-Hdr), 

where A^ is the only dimensionful parameter. To compute the A’s we adapt the NMSSMTools 
[60] code to the UMSSM. Here we do not consider pure UMSSM corrections from gauge 
contributions since the U(l)' gauge coupling is small compared to the Yukawa coupling 
of the top quark [15]. 
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The minimization conditions for the loop improved Higgs potential become 
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The corrected CP-even mass-squared matrix elements are 

(Xp.i = (xp,, + - 3A«| + As|) ^ 
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The corrected CP-odd mass-squared matrix elements read 
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which leads to the one-loop corrected pseudoscalar mass 
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where A567 = —A 5 sin 2/5 + AeC^ + Ayc^ and x = nsin2/5/2ns. 

Note that in NMSSMTools the CP-odd matrix M._ is defined in the basis of the two 
CP-odd bosons we must therefore perform the transformation T^Ai-T where 




-S/3 

0 ^ 

T = 


C/3 

0 


l 0 

0 

1/ 


(AS) 


as in [123]. 

Finally, the charged Higgs mass is corrected as 


{m%±f = (m^±)^ + 


A.^ + Ae^ + Ar^ 


'V2 


sin 2/3 


(A 4 -I- As)^. 


(A.7) 


References 

[1] CMS Collaboration, S. Chatrchyan et al, “Observation of a new boson at a mass 
of 125 GeV with the CMS experiment at the LHC,” Phys. Lett. B716 (2012) 
30-61, arXiv;1207.7235 [hep-ex] . 

[2] ATLAS Collaboration, G. Aad et al, “Observation of a new particle in the search 
for the Standard Model Higgs boson with the ATLAS detector at the LHC,” Phys. 
Lett. B716 (2012) 1-29, arXiv: 1207.7214 [hep-ex]. 

[3] CMS Collaboration, S. Chatrchyan et al, “Observation of a new boson with mass 
near 125 GeV in pp collisions at ^/s = 7 and 8 TeV,” JHEP 1306 (2013) 081, 
arXiv:1303.4571 [hep-ex]. 

[4] CMS Collaboration, V. Khachatryan et al., “Search for top-squark pairs decaying 
into Higgs or Z bosons in pp collisions at sqrt(s) = 8 TeV,” Phys. Lett. B736 
(2014) 371, arXiv:1405.3886 [hep-ex]. 

[5] CMS Collaboration, V. Khachatryan et al., “Searches for electroweak production 
of charginos, neutralinos, and sleptons decaying to leptons and W, Z, and Higgs 
bosons in pp collisions at 8 TeV,” Eur. Phys. J. C74 no. 9, (2014) 3036, 
arXiv:1405.7570 [hep-ex]. 

[6] ATLAS Collaboration, G. Aad et al., “Search for squarks and gluinos with the 
ATLAS detector in final states with jets and missing transverse momentum using 

= 8 TeV proton-proton collision data,” JHEP 1409 (2014) 176, 
arXiv:1405.7875 [hep-ex]. 

[7] ATLAS Collaboration, G. Aad et al., “Search for top squark pair production in 
final states with one isolated lepton, jets, and missing transverse momentum in 
^/s =8 TeV pp collisions with the ATLAS detector,” JHEP 1411 (2014) 118, 
arXiv:1407.0583 [hep-ex] . 

[8] L. J. Hall, D. Pinner, and J. T. Ruderman, “A Natural SUSY Higgs Near 126 
GeV,” JHEP 1204 (2012) 131, arXiv: 1112.2703 [hep-ph] . 


44 





[9] H. Baer, V. Barger, D. Mickelson, and M. Padeffke-Kirkland, “SUSY models 
under siege: LHC constraints and electroweak fine-tuning,” Phys. Rev. D89 
(2014) 115019, arXiv:1404.2277 [hep-ph] . 

[10] U. Ellwanger and C. Hugonie, “The semi-constrained NMSSM satisfying bounds 
from the LHC, LUX and Planck,” JHEP 1408 (2014) 046, arXiv: 1405.6647 
[hep-ph] . 

[11] A. Kaminska, G. G. Ross, K. Schmidt-Hoberg, and F. Staub, “A precision study 
of the fine tuning in the DiracNMSSM,” JHEP 1406 (2014) 153, 

arXiv:1401.1816 [hep-ph]. 

[12] M. Farina, M. Perelstein, and B. Shakya, “Higgs Couplings and Naturalness in 
A-SUSY,” JHEP 1404 (2014) 108, arXiv: 1310.0459 [hep-ph]. 

[13] P. Athron, M. Binjonaid, and S. F. King, “Fine Tuning in the Constrained 
Exceptional Supersymmetric Standard Model,” Phys. Rev. D87 no. 11, (2013) 
115023, arXiv:1302.5291 [hep-ph]. 

[14] M. Cvetic, D. A. Demir, J. Espinosa, L. Everett, and P. Langacker, “Electroweak 
breaking and the mu problem in supergravity models with an additional U(l),” 
Phys. Rev. D56 (1997) 2861, arXiv:hep-ph/9703317 [hep-ph]. 

[15] V. Barger, P. Langacker, H.-S. Lee, and G. Shaughnessy, “Higgs Sector in 
Extensions of the MSSM,” Phys. Rev. D73 (2006) 115010, 

arXiv:hep-ph/0603247 [hep-ph]. 

[16] ATLAS Collaboration, G. Aad et ai, “Search for high-mass dilepton resonances 
in pp collisions at ^/s = 8TeV with the ATLAS detector,” Phys. Rev. D90 no. 5, 
(2014) 052005, arXiv:1405.4123 [hep-ex]. 

[17] CMS Collaboration, V. Khachatryan et ai, “Search for physics beyond the 

standard model in dilepton mass spectra in proton-proton collisions at ^/s = 8 
TeV,” 1504 (2015) 025, arXiv: 1412.6302 [hep-ex]. 

[18] M. Frank and S. Mondal, “Light neutralino dark matter in U(l)' models,” Phys. 
Rev. D90 no. 7, (2014) 075013, arXiv: 1408.2223 [hep-ph]. 

[19] ATLAS Collaboration, “Measurements of the Higgs boson production and decay 
rates and coupling strengths using pp collision data at ^/s = 7 and 8 TeV in the 
ATLAS experiment,”, http://cds.cern.ch/record/2002212. 
ATLAS-CONF-2015-007, ATLAS-COM-CONF-2015-011. 

[20] CMS Collaboration, V. Khachatryan et ai, “Precise determination of the mass of 
the Higgs boson and tests of compatibility of its couplings with the standard 
model predictions using proton collisions at 7 and 8 TeV,” Eur. Phys. J. C75 

no. 5,(2015) 212, arXiv:1412.8662 [hep-ex]. 

[21] P. Athron, D. Harries, and A. Williams, “Z’ mass limits and the naturalness of 
supersymmetry,” arXiv: 1503.08929 [hep-ph] . 


45 


[22] M. Cvetic and P. Langacker, “Implications of Abelian extended gauge structures 
from string models,” Phys. Rev. D54 (1996) 3570-3579, arXiv:hep-ph/9511378 
[hep-ph] . 

[23] M. Cvetic and P. Langacker, “New gauge bosons from string models,” Mod. Phys. 
Lett. All (1996) 1247-1262, arXiv;hep-ph/9602424 [hep-ph], 

[24] G. Cleaver, M. Cvetic, J. R. Espinosa, L. L. Everett, and P. Langacker, 
“Classification of fiat directions in perturbative heterotic superstring vacua with 
anomalous U(l),” Nucl. Phys. B525 (1998) 3-26, arXiv;hep-th/9711178 
[hep-th] . 

[25] G. Cleaver, M. Cvetic, J. Espinosa, L. Everett, P. Langacker, et ai, “Physics 
implications of fiat directions in free fermionic superstring models 1. Mass 
spectrum and couplings,” Phys. Rev. D59 (1999) 055005, arXiv:hep-ph/9807479 
[hep-ph] . 

[26] G. Cleaver, M. Cvetic, J. Espinosa, L. Everett, P. Langacker, et ai, “Physics 
implications of fiat directions in free fermionic superstring models. 2. 
Renormalization group analysis,” Phys. Rev. D59 (1999) 115003, 

arXiv;hep-ph/9811355 [hep-ph], 

[27] P. Langacker, “Grand Unified Theories and Proton Decay,” Phys. Rept. 72 (1981) 
185. 

[28] D. London and J. L. Rosner, “Extra Gauge Bosons in E(6),” Phys. Rev. D34 
(1986) 1530. 

[29] P. Athron, S. King, D. Miller, S. Moretti, and R. Nevzorov, “Constrained 
Exceptional Supersymmetric Standard Model with a Higgs Near 125 GeV,” Phys. 
Rev. D86 (2012) 095003, arXiv: 1206.5028 [hep-ph], 

[30] K. Cheung, C.-T. Lu, and T.-C. Yuan, “Diphoton Rate of the 
Standard-Model-Like Higgs Boson in the Extra U(l) Extended MSSM,” Phys. 

Rev. D87 (2013) 075001, arXiv: 1212.1288 [hep-ph], 

[31] L. Basso and F. Staub, “Enhancing h 77 with staus in SUSY models with 
extended gauge sector,” Phys. Rev. D87 no. 1, (2013) 015011, arXiv: 1210.7946 
[hep-ph] . 

[32] M. Frank, L. Selbuz, L. Solmaz, and 1. Turan, “Higgs bosons in supersymmetric 
U(l)’ models with CP violation,” Phys. Rev. D87 no. 7, (2013) 075007, 
arXiv:1302.3427 [hep-ph] . 

[33] P. Athron, M. Mhlleitner, R. Nevzorov, and A. Williams, “Non-Standard Higgs 
Decays in U(l) Extensions of the MSSM,” JHEP 1501 (2015) 153, 

arXiv:1410.6288 [hep-ph], 

[34] G. Belanger, J. Da Silva, and A. Pukhov, “The Right-handed sneutrino as thermal 
dark matter in U(l) extensions of the MSSM,” JCAP 1112 (2011) 014, 

arXiv:1110.2414 [hep-ph]. 


46 


[35] L. Basso, B. O’Leary, W. Porod, and F. Staub, “Dark matter scenarios in the 
minimal SUSY B-L model,” JHEP 1209 (2012) 054, arXiv: 1207.0507 [hep-ph] . 

[36] M. Hirsch, W. Porod, L. Reichert, and F. Staub, “Phenomenology of the minimal 
supersymmetric U{1)b-l x U{1)r extension of the standard model,” Phys. Rev. 
D86 (2012) 093018, arXiv; 1206 .3516 [hep-ph]. 

[37] J. Camargo-Molina, B. O’Leary, W. Porod, and F. Staub, “The Stability Of 
R-Parity In Supersymmetric Models Extended By U{1)b-l” Phys. Rev. D88 
(2013) 015033, arXiv;1212.4146 [hep-ph]. 

[38] S. Di Chiara, V. Keus, and O. Lebedev, “Stabilizing the Higgs potential with a 
Z',” Phys. Lett. B744 (2015) 59-66, arXiv: 1412.7036 [hep-ph]. 

[39] S. Kraml, S. Kulkarni, U. Laa, A. Lessa, W. Magerl, et ai, “SModelS: a tool for 
interpreting simplified-model results from the LHC and its application to 
supersymmetry,” Eur. Phys. J. C74 (2014) 2868, arXiv; 1312.4175 [hep-ph]. 

[40] S. Kraml, S. Kulkarni, U. Laa, A. Lessa, V. Magerl, W. Magerl, 

D. Proschofsky-Spindler, M. Traub, and W. Waltenberger, “SModelS vl.O: a short 
user guide,” arXiv: 1412.1745 [hep-ph] . 

[41] M. A. Ajaib, B. Dutta, T. Ghosh, 1. Gogoladze, and Q. Shafi, “Neutralinos and 
Sleptons at the LHC in Light of Muon [g — 2)^,” arXiv; 1505.05896 [hep-ph] . 

[42] J. Kalinowski, S. King, and J. Roberts, “Neutralino Dark Matter in the USSM,” 
Ji/EP 0901 (2009) 066, arXiv: 0811.2204 [hep-ph]. 

[43] P. Langacker and J. Wang, “U(l)' symmetry breaking in supersymmetric E(6) 
models,” Phys. Rev. D58 (1998) 115010, arXiv:hep-ph/9804428 [hep-ph]. 

[44] V. Barger, P. Langacker, 1. Lewis, M. McCaskey, G. Shaughnessy, et ai, “Recoil 
Detection of the Lightest Neutralino in MSSM Singlet Extensions,” Phys. Rev. 
D75 (2007) 115002, arXiv:hep-ph/0702036 [HEP-PH]. 

[45] P. Langacker, “The Physics of Heavy Z' Gauge Bosons,” Rev. Mod. Phys. 81 
(2009) 1199-1228, arXiv ;0801. 1345 [hep-ph]. 

[46] B. O’Leary, W. Porod, and E. Staub, “Mass spectrum of the minimal SUSY B-L 
model,” JHEP 1205 (2012) 042, arXiv; 1112.4600 [hep-ph]. 

[47] M. E. Krauss, B. O’Leary, W. Porod, and E. Staub, “Implications of gauge kinetic 
mixing on Z’ and slepton production at the LHG,” Phys. Rev. D86 (2012) 055017, 
arXiv;1206.3513 [hep-ph]. 

[48] A. Leike, S. Riemann, and T. Riemann, “Z Z' mixing in presence of standard 
weak loop corrections,” arXiv;hep-ph/9808374 [hep-ph]. 

[49] J. Erler, P. Langacker, S. Munir, and E. Rojas, “Improved Constraints on Z' 
Bosons from Electroweak Precision Data,” JHEP 0908 (2009) 017, 

arXiv;0906.2435 [hep-ph] . 


47 


[50] D. Cerdeno, C. Hugonie, D. Lopez-Fogliani, C. Munoz, and A. Teixeira, 
“Theoretical predictions for the direct detection of neutralino dark matter in the 
NMSSM,” JHEP 0412 (2004) 048, arXiv:hep-ph/0408102 [hep-ph] . 

[51] Muon g-2 Collaboration, G. Bennett et al, “Final Report of the Muon E821 
Anomalous Magnetic Moment Measurement at BNL,” Phys. Rev. D73 (2006) 
072003, arXiv:hep-ex/0602035 [hep-ex], 

[52] B. L. Roberts, “Status of the Fermilab Muon (g — 2) Experiment,” Chin. Phys. 
C34 (2010) 741-744, arXiv; 1001.2898 [hep-ex], 

[53] S, Hesselbach, E, Franke, and H, Fraas, “Neutralinos in E(6) inspired 
supersymmetric U(l)’ models,” Eur. Phys. J. C23 (2002) 149-162, 

arXiv:hep-ph/0107080 [hep-ph], 

[54] S, Choi, H, Haber, J, Kalinowski, and P, Zerwas, “The Neutralino sector in the 
U(l)-extended supersymmetric standard model,” Nucl. Phys. B778 (2007) 

85-128, arXiv:hep-ph/0612218 [hep-ph], 

[55] B, de Carlos and J, Espinosa, “Cold dark matter candidate in a class of 
supersymmetric models with an extra U(l),” Phys. Lett. B407 (1997) 12-21, 

arXiv:hep-ph/9705315 [hep-ph], 

[56] V, Barger, C, Kao, P, Langacker, and H,-S, Lee, “Neutralino relic density in a 
supersymmetric 77(1)' model,” Phys. Lett. B600 (2004) 104-115, 

arXiv:hep-ph/0408120 [hep-ph], 

[57] S, King, S, Moretti, and R, Nevzorov, “Theory and phenomenology of an 
exceptional supersymmetric standard model,” Phys. Rev. D73 (2006) 035009, 

arXiv:hep-ph/0510419 [hep-ph], 

[58] P, Bechtle, O, Brein, S, Heinemeyer, O, Stal, T, Stefaniak, et al, 

“HiggsBounds — 4: Improved Tests of Extended Higgs Sectors against Exclusion 
Bounds from LEP, the Tevatron and the LHC,” Eur. Phys. J. C74 (2014) 2693, 
arXiv;1311,0055 [hep-ph], 

[59] P, Bechtle, S, Heinemeyer, O, Stal, T, Stefaniak, and G, Weiglein, ^'HiggsSignals 
Confronting arbitrary Higgs sectors with measurements at the Tevatron and the 
LHC,” Eur. Phys. J. C74 (2014) 2711, arXiv: 1305,1933 [hep-ph], 

[60] U, Ellwanger and C, Hugonie, “NMHDECAY 2,0: An Updated program for 
sparticle masses, Higgs masses, couplings and decay widths in the NMSSM,” 
Comput. Phys. Commun. 175 (2006) 290-303, arXiv:hep-ph/0508022 [hep-ph] 

[61] T, Sjostrand, S, Mrenna, and P, Z, Skands, “A Brief Introduction to PYTHIA 
8,1,” Comput. Phys. Commun. 178 (2008) 852-867, arXiv:0710,3820 [hep-ph], 

[62] A, Martin, W, Stirling, R, Thorne, and G, Watt, “Parton distributions for the 
LHG,” Eur. Phys. J. C63 (2009) 189-285, arXiv:0901,0002 [hep-ph]. 


48 


[63] E. Accomando, A. Belyaev, J. Fiaschi, K. Mimasu, S. Moretti, et ai, 
“Forward-Backward Asymmetry as a Discovery Tool for Z’ Bosons at the LHC,” 
arXiv:1503.02672 [hep-ph] . 

[64] J. Gunion, L. Roszkowski, and H. Haber, “Z' mass limits, masses and coupling of 
higgs bosons, and Z' decays in an Eq superstring based model,” Phys. Lett. B189 
(1987) 409. 

[65] T. Gherghetta, T. A. Kaeding, and G. L. Kane, “Supersymmetric contributions to 
the decay of an extra Z boson,” Phys. Rev. D57 (1998) 3178-3181, 

arXiv:hep-ph/9701343 [hep-ph], 

[66] C.-F. Chang, K. Cheung, and T.-C. Yuan, “Supersymmetric Decays of the Z' 
Boson,” JHEP 1109 (2011) 058, arXiv; 1107.1133 [hep-ph], 

[67] G. Corcella, “Phenomenology of supersymmetric Z’ decays at the Large Hadron 
Collider,” arXiv: 1412.6831 [hep-ph] . 

[68] K. Babu, C. F. Kolda, and J. March-Russell, “Leptophobic U(l)’s and the R(6) - 
R(c) crisis,” Phys. Rev. D54 (1996) 4635-4647, arXiv;hep-ph/9603212 
[hep-ph] . 

[69] Particle Data Group Collaboration, K. Olive et ai, “Review of Particle 
Physics,” Chin. Phys. C38 (2014) 090001. 

[70] G. Belanger, M. Heikinheimo, and V. Sanz, “Model-Independent Bounds on 
Squarks from Monophoton Searches,” JHEP 1208 (2012) 151, arXiv: 1205.1463 
[hep-ph] . 

[71] ATLAS Collaboration, G. Aad et ai, “Search for new phenomena in events with 
a photon and missing transverse momentum in pp collisions at ^/s = 8 TeV with 
the ATLAS detector,” Phys. Rev. D91 no. 1, (2015) 012008, arXiv: 1411.1559 
[hep-ex] . 

[72] H. K. Dreiner, M. Kramer, and J. Tattersall, “How low can SUSY go? Matching, 
monojets and compressed spectra,” Europhys. Lett. 99 (2012) 61001, 

arXiv:1207.1613 [hep-ph], 

[73] G. Belanger, F. Boudjema, A. Pukhov, and A. Semenov, “micrOMEGAs4.1: two 
dark matter candidates,” Comput. Phys. Commun. 192 (2015) 322-329, 

arXiv:1407.6129 [hep-ph], 

[74] W. Beenakker, R. Hopker, M. Spira, and P. Zerwas, “Squark and gluino 
production at hadron colliders,” Nucl. Phys. B492 (1997) 51-103, 
arXiv;hep-ph/9610490 [hep-ph], 

[75] W. Beenakker, M. Kramer, T. Plehn, M. Spira, and P. Zerwas, “Stop production 
at hadron colliders,” Nucl. Phys. B515 (1998) 3-14, arXiv:hep-ph/9710451 
[hep-ph] . 


49 


[76] A. Kulesza and L. Motyka, “Threshold resummation for squark-antisquark and 
gluino-pair production at the LHC,” Phys. Rev. Lett. 102 (2009) 111802, 
arXiv:0807.2405 [hep-ph] . 

[77] A. Kulesza and L. Motyka, “Soft gluon resummation for the production of 
gluino-gluino and squark-antisquark pairs at the LHC,” Phys. Rev. D80 (2009) 
095004, arXiv;0905.4749 [hep-ph] . 

[78] W. Beenakker, S. Brensing, M. Kramer, A. Kulesza, E. Laenen, et al, “Soft-gluon 
resummation for squark and gluino hadroproduction,” JHEP 0912 (2009) 041, 
arXiv;0909.4418 [hep-ph], 

[79] W. Beenakker, S. Brensing, M. Kramer, A. Kulesza, E. Laenen, et al, 
“Supersymmetric top and bottom squark production at hadron colliders,” JHEP 
1008 (2010) 098, arXiv; 1006.4771 [hep-ph]. 

[80] W. Beenakker, S. Brensing, M. Kramer, A. Kulesza, E. Laenen, et al, “Squark 
and Gluino Hadroproduction,” Int. J. Mod. Phys. A26 (2011) 2637-2664, 
arXiv;1105.1110 [hep-ph]. 

[81] DO Collaboration, V. M. Abazov et al, “Search for charged massive long-lived 
particles at ^/s = 1.96 TeV,” Phys. Rev. D87 no. 5, (2013) 052011, 

arXiv;1211.2466 [hep-ex]. 

[82] ATLAS Collaboration, G. Aad et al, “Searches for heavy long-lived charged 
particles with the ATLAS detector in proton-proton collisions at ^/s = 8 TeV,” 
JHEP 1501 (2015) 068, arXiv; 1411.6795 [hep-ex]. 

[83] CMS Collaboration, S. Chatrchyan et al, “Searches for long-lived charged 
particles in pp collisions at ■s/s=l and 8 TeV,” JHEP 1307 (2013) 122, 
arXiv;1305.0491 [hep-ex]. 

[84] CMS Collaboration, V. Khachatryan et al., “Constraints on the pMSSM, AMSB 
model and on other models from the search for long-lived charged particles in 
proton-proton collisions at -y/s = 8 TeV,” arXiv; 1502.02522 [hep-ex]. 

[85] LHCb Collaboration, R. Aaij et al., “Implications of LHCb measurements and 
future prospects,” Eur. Phys. J. C73 no. 4, (2013) 2373, arXiv; 1208.3355 
[hep-ex] . 

[86] N. Bernal, M. Losada, and E. Mahmoudi, “Flavour physics constraints in the 
BMSSM,” JHEP 1107 (2011) 074, arXiv; 1104.5395 [hep-ph]. 

[87] A. Arbey, M. Battaglia, F. Mahmoudi, and D. M. Santos, “Supersymmetry 
confronts Bg —?• iP~ii~'. Present and future status,” Phys. Rev. D87 (2013) 035026, 

arXiv;1212.4887 [hep-ph]. 

[88] F. Domingo and U. Ellwanger, “Updated Constraints from B Physics on the 
MSSM and the NMSSM,” JHEP 0712 (2007) 090, arXiv; 0710.3714 [hep-ph]. 

[89] J. Da Silva, “Supersymmetric Dark Matter candidates in light of constraints from 
collider and astroparticle observables,” arXiv; 1312.0257 [hep-ph] . 


50 


[90] Particle Data Group Collaboration, J. Beringer et ai, “Review of Particle 
Physics (RPP),” Phys. Rev. D86 (2012) 010001. 

[91] Heavy Flavor Averaging Group, www.slac.stanford.edu/xorg/hfag/rare/ 
2013/radll/0UTPUT/HTML/radll_table7.html. 

[92] Heavy Flavor Averaging Group. 
www.slac.stanford.edu/xorg/hfag/rare/2013/radll/btosg.pdf. 

[93] CMS, LHCb Collaboration, V. Khachatryan et aL, “Observation of the rare 

> /r+/r“ decay from the combined analysis of CMS and LHCb data,” Nature 
(2015) , arXiv:1411.4413 [hep-ex]. 

[94] Heavy Flavor Averaging Group. 

WWW .slac.Stanford.edu/xorg/hf ag/osc/PDG_2014/#DMS. 

[95] Heavy Flavor Averaging Group. 

WWW .slac.Stanford.edu/xorg/hf ag/osc/PDG_2014/#DMD. 

[96] T. Aoyama, M. Hayakawa, T. Kinoshita, and M. Nio, “Complete Tenth-Order 
QED Contribution to the Muon g-2,” Phys. Rev. Lett. 109 (2012) 111808, 
arXiv;1205.5370 [hep-ph] . 

[97] Planck Collaboration, P. Ade et a/., “Planck 2015 results. XHI. Cosmological 
parameters,” arXiv: 1502.01589 [astro-ph. CO] . 

[98] L. J. Hall, K. Jedamzik, J. March-Russell, and S. M. West, “Freeze-In Production 
of FIMP Dark Matter,” JHEP 1003 (2010) 080, arXiv:0911.1120 [hep-ph] . 

[99] X. Chu, T. Hambye, and M. H. Tytgat, “The Four Basic Ways of Creating Dark 
Matter Through a Portal,” JCAP 1205 (2012) 034, arXiv: 1112.0493 [hep-ph]. 

[100] F. Franke and S. Hesselbach, “Production of singlino dominated neutralinos in 
extended supersymmetric models,” Phys. Lett. B526 (2002) 370-378, 

arXiv:hep-ph/0111285 [hep-ph]. 

[101] D. Suematsu, “Singlino dominated LSP as CDM candidate in supersymmetric 
models with an extra U(l),” Phys. Rev. D73 (2006) 035010, 

arXiv:hep-ph/0511299 [hep-ph]. 

[102] S. Nakamura and D. Suematsu, “Supersymmetric extra U(l) models with a 
singlino dominated LSP,” Phys. Rev. D75 (2007) 055004, arXiv:hep-ph/0609061 
[hep-ph] . 

[103] LUX Collaboration Collaboration, D. Akerib et al, “First results from the LUX 
dark matter experiment at the Sanford Underground Research Facility,” Phys. 

Rev. Lett. 112 no. 9, (2014) 091303, arXiv: 1310.8214 [astro-ph.CO] . 

[104] Fermi-LAT Collaboration, M. Ackermann et ai, “Searching for Dark Matter 
Annihilation from Milky Way Dwarf Spheroidal Galaxies with Six Years of 
Fermi-LAT Data,” arXiv: 1503.02641 [astro-ph. HE] . 


51 


[105] ATLAS, CDF, CMS, DO Collaboration, “First combination of Tevatron and 
LHC measurements of the top-quark mass,” arXiv; 1403.4427 [hep-ex] . 

[106] CMS Collaboration, “Search for New Physics in the Multijets and Missing 
Momentum Final State in Proton-Proton Collisions at 8 TeV,”. 

http: //cds. cern. ch/record/1563156. CMS-PAS-SUS-13-012. 

[107] CMS Collaboration, “Search for supersymmetry in hadronic final states using 
Mt 2 with the CMS detector at sqrt(s) = 8 TeV,”. 

http: //cds. cern. ch/record/1646394. CMS-PAS-SUS-13-019. 

[108] ATLAS Collaboration, “Search for strongly produced superpartners in final states 
with two same sign leptons with the ATLAS detector using 21 fb“^ of 
proton-proton collisions at ^/s = 8 TeV,”. http;//cds.cern.ch/record/2002212. 
ATLAS-CONF-2013-007, ATLAS-COM-CONF-2013-006. 

[109] ATLAS Collaboration, “Search for strong production of supersymmetric particles 
in final states with missing transverse momentum and at least three b-jets using 
20.1 fb“^ of pp collisions at ^/s = 8 TeV with the ATLAS Detector,”. 

https://cds.cern.ch/record/1557778. ATLAS-CONF-2013-061, 
ATLAS-COM-CONF-2013-071. 

[110] CMS Collaboration, “Search for supersymmetry in final states with missing 
transverse energy and 0, 1, 2, 3, or at least 4 b-quark jets in 8 TeV pp collisions 
using the variable «r,”. http://cds.cern.ch/record/1494483. 
CMS-PAS-SUS-12-028. 

[111] CMS Collaboration, S. Chatrchyan et ai, “Search for supersymmetry in pp 
collisions at yi=8 TeV in events with a single lepton, large jet multiplicity, and 
multiple b jets,” Phys. Lett. B733 (2014) 328-353, arXiv; 1311.4937 [hep-ex]. 

[112] CMS Collaboration, “Exclusion limits on gluino and top-squark pair production 
in natural SUSY scenarios with inclusive razor and exclusive single-lepton searches 
at 8 TeV,”. http;//cds.cern.ch/record/1745586. CMS-PAS-SUS-14-011. 

[113] CMS Collaboration, S. Chatrchyan et ai, “Search for gluino mediated bottom- 
and top-squark production in multijet final states in pp collisions at 8 TeV,” Phys. 
Lett. B725 (2013) 243-270, arXiv: 1305.2390 [hep-ex]. 

[114] ATLAS Collaboration. https://atlas.web.cern.ch/Atlas/GROUPS/PHYSICS/ 
C0NFN0TES/ATLAS-C0NF-2013-061/fig_12c.png. 

[115] CMS Collaboration, V. Khachatryan et ai, “Search for supersymmetry using 
razor variables in events with 5-tagged jets in pp collisions at ^/s = 8 TeV,” Phys. 
Rev. D91 (2015) 052018, arXiv: 1502.00300 [hep-ex]. 

[116] CMS Collaboration, “Search for direct production of bottom squark pairs,”, 
http: //cds. cern. ch/record/1693164. CMS-PAS-SUS-13-018. 

[117] CMS Collaboration, S. Chatrchyan et ai, “Search for top-squark pair production 
in the single-lepton final state in pp collisions at ^/s = 8 TeV,” Eur. Phys. J. C73 
no. 12, (2013) 2677, arXiv:1308.1586 [hep-ex]. 


52 


[118] ATLAS Collaboration, G. Aad et ai, “Search for direct third-generation squark 
pair production in final states with missing transverse momentum and two 6-jets 
in -y/s = 8 TeV pp collisions with the ATLAS detector,” JHEP 1310 (2013) 189, 
arXiv:1308.2631 [hep-ex], 

[119] L. Edelhuser, J. Heisig, M. Krmer, L. Oymanns, and J. Sonneveld, “Constraining 
supersymmetry at the LHC with simplified models for squark production,” JHEP 
1412 (2014) 022, arXiv:1410.0965 [hep-ph] . 

[120] S. Kraml et ai, “Smodels vl.l,”. in preparation. 

[121] ATLAS Collaboration, G. Aad et ai, “Search for direct production of charginos, 
neutralinos and sleptons in final states with two leptons and missing transverse 
momentum in pp collisions at \/s = 8 TeV with the ATLAS detector,” JHEP 
1405 (2014) 071, arXiv; 1403.5294 [hep-ex]. 

[122] C. Arina, M. E. C. Catalan, S. Kraml, S. Kulkarni, and U. Laa, “Constraints on 
sneutrino dark matter from LHC Run 1,” arXiv: 1503.02960 [hep-ph] . 

[123] G. Belanger, E. Boudjema, C. Hugonie, A. Bukhov, and A. Semenov, “Relic 
density of dark matter in the NMSSM,” JCAP 0509 (2005) 001, 

arXiv:hep-ph/0505142 [hep-ph]. 

[124] G. Belanger, D. Ghosh, R. Godbole, and S. Kulkarni, “Light stop in the MSSM 
after LHC Run 1,” arXiv: 1506.00665 [hep-ph]. 

[125] CMS Collaboration, S. Chatrchyan et ai, “Search for invisible decays of Higgs 
bosons in the vector boson fusion and associated ZH production modes,” Eur. 
Phys. J. C74 (2014) 2980, arXiv: 1404.1344 [hep-ex]. 

[126] ATLAS Collaboration, “Search for an Invisibly Decaying Higgs Boson Produced 
via Vector Boson Fusion in pp Collisions at ^/s = 8 TeV using the ATLAS 
Detector at the LHC,”. http;//cds.cern.ch/record/2002121. 
ATLAS-CONF-2015-004, ATLAS-COM-CONF-2015-004. 

[127] J. Bernon, B. Dumont, and S. Kraml, “Status of Higgs couplings after run 1 of the 
LHC,” Phys. Rev. D90 (2014) 071301, arXiv: 1409.1588 [hep-ph]. 

[128] D. Ghosh, R. Godbole, M. Guchait, K. Mohan, and D. Sengupta, “Looking for an 
Invisible Higgs Signal at the LHC,” Phys. Lett. B725 (2013) 344-351, 

arXiv:1211.7015 [hep-ph]. 

[129] A. Djouadi, L. Maiani, A. Polosa, J. Quevillon, and V. Riquer, “Fully covering the 
MSSM Higgs sector at the LHC,” arXiv: 1502.05653 [hep-ph] . 

[130] M. Bisset, J. Li, N. Kersting, R. Lu, F. Moortgat, et al, “Four-lepton LHC events 
from MSSM Higgs boson decays into neutralino and chargino pairs,” JHEP 0908 
(2009) 037, arXiv:0709.1029 [hep-ph]. 

[131] A. J. Williams, C. Boehm, S. M. West, and D. A. Vasquez, “Regenerating WIMPs 
in the Light of Direct and Indirect Detection,” Phys. Rev. D86 (2012) 055018, 
arXiv:1204.3727 [hep-ph] . 


53 


[132] T. Cohen, M. Lisanti, A. Pierce, and T. R. Slatyer, “Wino Dark Matter Under 
Siege,” JCAP 1310 (2013) 061, arXiv: 1307.4082. 

[133] G. Belanger, C. Boehm, M. Cirelli, J. Da Silva, and A. Pukhov, “PAMELA and 
FERMI-LAT limits on the neutralino-chargino mass degeneracy,” JCAP 1211 
(2012) 028, arXiv:1208.5009 [hep-ph] . 

[134] M. L. Graesser and J. Shelton, “Hunting Mixed Top Squark Decays,” Phys. Rev. 
Lett. Ill no. 12, (2013) 121802, arXiv; 1212.4495 [hep-ph]. 


54 


